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ABSTRACT 


Flows  in  ducts  with  curvature  manifest  secondary  motions  driven  by 
lateral  curvature  of  the  main  flow.  In  addition,  cross-stream  motions  can 
arise  in  non-circular  ducts  due  to  gradients  of  the  Reynolds  stresses  in  the 
plane  perpendicular  to  the  main  flow.  Observations  from  this  study  suggest 
that  in  a  strongly  curved  180°  bend  of  square  cross-section,  the  latter 
mechanism  contributes  negligibly  to  the  secondary  motion.  New,  non-intru- 
slve  measurements  of  the  mean  flow  in  this  configuration  are  in  qualitative 
agreement  with  earlier  results  from  inviscid  flow  analysis.  Also,  measure¬ 
ments  of  the  turbulent  stresses  display  previously  undocumented  anisotropic 
characteristics  arising  from  shearing  motions  induced  in  the  core  of  the 
flow. 

Laminar  and  turbulent  flow  predictions  have  been  made  for  the  experi¬ 
mental  configuration.  Unlike  earlier  work,  turbulence  model  deficiencies 
are  more  clearly  revealed  by  the  higher  order  accuracy  of  the  finite  dif¬ 
ferencing  scheme  used  and  the  implementation  of  a  partially-parabollc  calcu¬ 
lation  algorithm.  It  is  found  that,  with  the  problem  of  numerical  diffusion 
relieved,  false  physical  diffusion  in  the  turbulence  model  is  the  main  cause 
for  the  differences  observed  between  measurements  and  computations.  Solutions 
of  the  turbulent  flow  field  in  the  bend  and  downstream  tangent  show  an  insensi 
tlvlty  to  variations  in  the  specification  of  the  inlet  plane  boundary  condi¬ 
tion  for  streamwise  vortlcity.  An  explanation  for  this  finding  is  given  In 
terms  of  an  approximate  equation  for  the  evolution  of  streamwise  vortlcity. 

Knowledge  of  the  type  of  flow  investigated  here  is  essential  for  testing 
and  advancing  new  turbulence  modeling  concepts.  As  examples,  in  addition  to 
the  main  secondary  flow,  laminar  flow  predictions  reveal  a  Dean-type  vortex 
structure  at  the  concave  wall  in  the  bend;  a  second  vortex  on  the  convex  wall; 


and  a  third  corner  vortex  on  the  convex  wall-side  of  the  bend.  The  third 
vortex  is  associated  with  the  phenomenon)  of  wall  stall  in  cascade  flows. 

In  the  downstream  tangent,  turbulent  flow  measurements  show  drastic 
reductions  of  the  cross-stream  motion  taking  place  in  less  than  5  hydraulic 
diameters,  but  the  main  flow  recovers  only  slowly  from  the  effects  of  the 
bend.  Emphasis  should  be  placed  on  the  development  of  higher  order  model 
closures  for  predicting  the  type  of  highly  anisotropic  three-dimensional 
turbulent  flow  examined  here. 
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local  radius  of  curvature  of  a  streamline 

r.  +  r 
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transverse  or  radial  coordinate  in  bend 

r1 
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(=  (r  -  r^)/(r0  ■  <^))  normalized  transverse  or  radial  coordinate 

streamwise  or  longitudinal  coordinate  in  bend 
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1 .  INTRODUCTION 

1.1  Secondary  Motions  in  Ducted  Flows 

Flows  in  ducts  with  curvature  experience  pressure  induced  cross-stream 
(secondary)  motions  of  academic  interest,  and  having  very  significant  conse¬ 
quences  in  many  industrial  applications  involving  heat  and  mass  transfer. 
Further  to  the  motions  induced  by  the  transverse  pressure  gradient,  gradients 
of  the  Reynolds  stresses  in  the  cross-stream  plane  can  also  produce  a  secondary 
flow.  The  latter  flow  is  sometimes  referred  to  as  secondary  flow  of  the 
"second  kind"  and,  unlike  the  pressure-gradient  driven  flow  of  the  "first 
kind",  does  not  arise  in  the  laminar  flow  regime.  Typically,  motions  of 
the  first  kind  are  10-30  times  faster  than  motions  of  the  second  kind, 
suggesting  that  the  cross-stream  flow  in  a  curved  duct  Is  principally 
dominated  by  mean  pressure  forces.  The  theory  of  secondary  flow  induced 
by  lateral  curvature  of  the  main  flow  and  by  gradients  of  the  Reynolds 
stresses  respectively  has  been  reviewed  recently  by  Johnston  [1]. 

Experimental  measurements  of  develop  ng  mean  flow  and  turbulence  char¬ 
acteristics  in  a  90°  bend  of  square  cross-section  have  been  reported  by 
Humphrey,  Whitelaw  and  Yee  [2]  and  by  Taylor,  Whltelaw  and  Ylanneskis  [3]. 

These  studies  were  performed  using  non-disturbing  laser-Doppler  velocimetry 
techniques.  In  [3]  thin  boundary  layers  were  purposefully  induced  at  the 
bend  Inlet  plane,  whereas  In  [2]  essentially  fully  developed  straight  duct 
flow  was  presented  as  the  Inlet  condition.  The  most  marked  differences 
between  the  results  of  these  two  studies  were,  In  general,  the:  a)  stronger 
secondary  motions,  and  b)  higher  levels  of  turbulence,  experienced  by  the 
flow  with  thicker  inlet  boundary  layers.  In  a  study  also  addressing  the 
influence  of  inlet  plane  boundary  layer  thickness  on  secondary  motions,  Bruun 
[4]  has  made  Pitot  tube  measw  ~'*>ts  the  mean  velocity  and  pressure 
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distributions  in  a  120°  bend  of  rectangular  cross-section.  In  [4]  a  limited 
amount  of  total  pressure  field  data  is  reported  for  three  locations  in  a 
straight  duct  section  downstream  of  the  bend.  Detailed  measurements  of  the 
much  weaker  Reynolds  stress  driven  secondary  motions  for  the  case  of  developing 
flow  in  a  duct  of  square  cross-section  have  been  reported  by  Melling  and 
White! aw  [5],  and  in  a  series  of  papers  by  Gessner  et  al . ,  the  latest 
reference  being  [6].  It  is  noteworthy  that  the  flow  field  measurements  in 
[2,3]  do  not  show  the  symmetric  corner  vortex  structures  characteristic  of 
straight  duct  flow.  This  supports  the  notion  that  in  the  bend  sections  of 
these  studies  the  cross-stream  flow  was  influenced  more  stongly  by  the 
transverse  mean  pressure-gradient  than  by  gradients  of  the  Reynolds  stresses. 

Curved  duct  turbulent  flow  calculations  based  on  the  numerical  solution 
of  the  Navier-Stokes  equations  have  been  made  by  Humphrey,  Whitelaw  and  Yee 
[2]  for  their  experimental  configuration,  and  more  recently  by  Buggeln,  Briley 
and  McDonald  [7]  for  the  flow  measured  by  Taylor,  Whitelaw  and  Yianneskis  [3]. 

In  [2]  a  two-equation  k-c  model  of  turbulence  was  used  while  in  [7]  the  authors 
...iployed  the  one-equation  model  investigated  by  Shamroth  and  Gibeling  [8].  In 
general,  both  studies  show  good  agreement  between  measurements  and  calculations 
of  the  mean  velocity  components  up  to  a  bend  angle  of  approximately  45°,  but  as 
of  this  point  significant  deviations  appear.  For  bend  angles 
larger  than  45°  calculations  in  [7]  are  in  better  accord  with  the  measurements 
than  corresponding  results  in  [2],  This  is  attributed  primarily  to  the  thinner 
boundary  layer  inlet  condition  present  in  [3],  which  was  responsible  for  the 
production  of  weaker  secondary  motions  than  were  observed  and  predicted  in  [2]. 
It  should  be  remarked  that  cost  considerations  dictated  the  use  of  coarse  grids 
In  both  of  the  above  studies  and,  since  significant  levels  of  numerical  diffusion 
probably  affected  the  results,  there  is  an  inadequate  basis  for  judging  the 
merits  and  demerits  of  the  turbulence  models  used  respectively  by  the  authors. 


A  review  of  computational  methods  for  internal  flows  has  been  carried  out 
by  McNally  and  Sockol  [9].  The  following  is  a  brief  summary  relating  to 
numerical  techniques  of  relevance  to  this  work.  Within  the  class  of  numerical 
methods  which  neglects  streamwise  diffusion  is  the  partially-parabol ic  proce¬ 
dure  proposed  by  Pratap  and  Spalding  [10]  for  calculating  flow  in  curved  ducts. 

In  this  method,  iterated  forward  marching  sweeps  of  the  three-dimensional  flow 
field  are  performed,  yielding  solutions  of  the  momentum  equations  and  corrections 
to  the  pressure  field,  until  local  and  global  continuity  are  attained.  This 
type  of  numerical  procedure  was  used  by  Chang  et  al  [11]  to  predict  the  flow  of 
[2]  after  modifying  the  procedure  to  include  quadratic  upstream  interpolation 
of  convection  terms  as  in  Han  et  al  [12]  to  reduce  numerical  diffusion.  In 
general,  better  agreement  was  found  in  [11]  than  in  [2]  between  mean  velocity 
measurements  and  calculations  as  a  consequence  of  the  numerical  improvements 
to  the  calculation  scheme. 

The  advantages  of  partially-parabol ic  procedures  over  fully-elliptic  ones 
are  mainly  two.  First,  they  require  considerably  less  computer  memory  since, 
in  a  three-dimensional  flow  calculation,  pressure  is  the  only  variable  requiring 
full  three-dimensional  storage.  The  second  advantage  derives  from  the  first. 

Due  to  the  reduced  memory  requirements,  with  a  partial ly-parabolic  procedure 
it  becomes  feasible  to  compute  three-dimensional  flows  on  calculation  grids 
refined  to  the  point  where  streamwise  numerical  diffusion  is  a  substantially 
weaker  source  of  numerical  inaccuracy.  A  disadvantage  of  the  partially-para¬ 
bol  ic  procedure  is  the  requirement  of  calculation  times  comparable  to  those 
of  an  elliptic  procedure. 

Inviscid  flow  approximations  for  the  calculation  of  flows  in  curved 
ducts  have  been  proposed  by,  for  example,  Briley  and  McDonal  [13].  The 
numerical  procedure  developed  by  these  authors  represents  an  extension  of 
viscous  forward  marching  methods  in  that  it  accounts  approximately  for 
transverse  variations  in  streanwise  pressure  gradient.  While  this  approach 


appears  particularly  promising  fo r  fast  and  relatively  inexpensive  curved  flow 
computations,  work  continues  on  its  development  and  application;  see  Krekovsky, 
Briley  and  McDonald  [14]. 


In  contrast  to  the  pressure-dominated  flows  arising  in  curved  ducts,  the 
difficulties  associated  with  predicting  Reynolds  stress  driven  secondary 
motions  in  ducts  of  non-circular  cross-section  are  principally  related  to  the 
accurate  modeling  of  cross-stream  turbulent  flow  anisotropy.  Brundrett  and 
Baines  [15]  have  shown  that  the  turbulence  term  contributing  most  strongly  to 
mean  streamwise  vorticity  in  steady,  incompressible,  constant  property  flow  is 

«  — 7 

(ul7  -  u,  ).  In  a  later  study,  Perkins  [16]  also  demonstrated  the  impor- 
ojoZ  y  Z  99 

_ 

tance  of  the  term  i-2-* - *•)  u  u  .  In  these  terms,  y,  z,  u  and  u  are  the 

3z  3y  y  y 

cross-stream  coordinates  and  velocity  fluctuations  respectively.  Since 

straight  duct  flows  are  correctly  described  by  parabolic  (boundary 
layer)  equations,  they  do  not  pose  serious  impediments  to  accurate  numerical 
solutions;  the  limitations  are  related  to  the  turbulence  models  used.  It  can 
readily  be  shown  that  turbulence  models  based  on  the  notion  of  an  isotropic 
turbulent  viscosity,  with  the  implication  that  the  principal  axes  of  stress 
and  strain  are  co-aligned,  cannot  account  for  the  generation  of  stress-driven 
secondary  motions.  In  their  recent  study,  Naot  and  Rodi  [17]  have  reviewed 
briefly  various  investigations  employing  Reynolds  stress  level  closures  to  pre¬ 
dict  these  motions.  They  conclude  that  the  essential  basic  flow  features 
can  be  predicted  qualitatively  with  relatively  simple  algebraic  stress  closures 
It  is  noteworthy  that,  except  for  the  full  Reynolds  stress  closure  predictions 
by  Reece  [18],  all  numerical  studies  of  the  problem  have  neglected  convection 
and  diffusion  contributions  to  the  balance  of  u.u . ,  thus  implying  the  condition 

*  J 

of  an  equilibrium  turbulent  flow. 
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1 .2  The  Problem  of  Interest 

The  flow  configuration  of  interest  to  this  study  consists  of  a  180°  bend 
of  square  cross-section  preceded  and  followed  by  straight  duct  sections.  In¬ 
side  the  bend  section  the  flow  pattern  is  established  primarily  by  pressure- 
force  imbalances.  Within  the  downstream  tangent  the  flow  undergoes  a  relaxa¬ 
tion  driven  by  turbulence  diffusion  and  redistribution  processes.  The 
flow  configuration  is  directly  relevant  to  heat  exchangers,  and  is  typical  of 
industrial  components  used  for  coal-slurry  transport.  The  flow  patterns 
established  in  and  downstream  of  the  bend  will  strongly  influence  the  detailed 
heat  transfer  characteristics  dictating  the  overall  performance  of  heat 
exchangers.  Knowledge  of  the  flow  and  of  its  turbulence  characteristics  is 
also  necessary  for  understanding  mechanisms  of  particle  transport  and  surface 
wear  by  slurries. 

1 .3  Objectives  of  this  Study 

The  objectives  of  this  study  have  been  mainly  two: 

1.  To  obtain  detailed  measurements  by  means  of  laser-Doppler  velocimetry. 
of  the  mean  flow  and  turbulence  characteristics  in  the  curved  and  down¬ 
stream  sections  of  the  above  configuration. 

In  contrast  to  earlier  work  [2-6],  the  present  configuration  involves 
the  combined  effects  of  a  180°  bend  and  downstream  tangent.  These  are 
responsible  for  inducing  a  complex  flow  of  very  distinct  and  previously 
undocumented  velocity  characteristics.  The  body  of  experimental  data  offered 
here  should  be  of  value  for  extending  the  theoretical  base  required  for 
developing  and  testing  more  comprehensive  three-dimensional  turbulence 
model  formulations. 

2.  To  test  the  capability  of  the  two  equation  (k-r)  model  of  turbulence 
to  predict  the  mean  features  of  the  flow  configuration. 
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Elliptic  calculations  in  [2]  of  the  turbulent  flow  in  a  90°  bend 
of  square  cross-section  yielded  results  in  poor  agreement  with  the 
measurements  for  bend  angles  larger  than  45°.  However,  in  that  study 
it  was  not  possible  to  separate  satifactorily  the  turbulence  model 
failings  from  inaccuracies  produced  by  numerical  diffusion.  Subse¬ 
quent  work  by  Chang  et  al  [11]  using  a  partially-parabolic  scheme 
and  higher  order  differencing  of  convective  terms  yielded  substantially 
improved  predictions.  Their  numerical  method  is  employed  here  to 
establish  more  precisely  the  limitations  of  the  k-e  turbulence  model 
for  predicting  strongly  curved  and  relaxing  three-dimensional  confined 
flows. 

1 .4  Outline  of  the  Paper 

The  remainder  of  this  paper  is  divided  into  four  sections.  In  order, 
these  discuss:  a)  The  apparatus,  instrumentation,  experimental  methodology 
and  measurement  error;  b)  The  governing  equations,  turbulence  model,  flow 
field  boundary  conditions,  finite  differencing  practices,  solution  algorithm 
and  calculated  test  cases;  c)  The  measured  and  calculated  results.  Including 
laminar  flow  calculations;  d)  The  major  conclusions  of  this  work  with  some 
specific  recommendations. 
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2.  THE  EXPERIMENT 

2.1  Flow  System  and  Instrumentation 

The  experimental  system  was  composed  of:  a  water  rig,  of  which  the  most 
important  component  was  the  flow  test  section;  a  laser-Doppler  velocimeter  and 
its  associated  electronic  instrumentation;  and  a  PDP  11/34  Digital  Equipment 
Corporation  minicomputer.  A  description  of  the  experimental  apparatus  and 
instrumentation  has  already  been  given  in  [19].  A  summary  of  that 
description,  including  later  modifications  to  tne  apparatus,  is  provided  here. 

The  basic  components  of  the  flow  test  section  are  shown  schematically 
in  Figure  1,  and  comprised  two  straight  ducts  and  a  bend  of  square  cross-sections 
The  tangents  were  each  31  hydraulic  diameters  long  and  were  respectively  attached 
to  the  0°  (inlet)  and  180°  (outlet)  plane  of  the  bend.  The  ratio  of  bend  mean 
radius  of  curvature  to  hydraulic  diameter  was  RC/DH  a  3.35.  The  bend  component 
was  constructed  by  machining  an  open,  curved  channel  of  mean  radius  of  curva¬ 
ture  Rc  ■  14.92  cm  (±  0.02  cm)  into  one  of  the  faces  of  a  large,  flat,  solid 
piece  of  plexiglass  7.6  cm  thick.  A  plexiglass  plate  1.27  cm  thick  was  machined 
to  fit  snugly  as  a  lid  over  the  open  channel  piece,  so  yielding  an  enclosed 
curved  duct  shape  of  cross-section  dimensions  4.45  *  4.45  cm  (±  0.02  x  0.02  cm  ) 
An  0-ring  seal  placed  between  the  open  channel  piece  and  the  lid  prevented 
water  leaks.  The  modular  form  of  this  construction  made  disassembly  easy  for 
surface  cleaning  purposes.  The  upstream  and  downstream  tangents  were  each 
138  cm  (±  0.02  cm)  long  and  were  also  constructed  from  plexiglass,  with  flat 
walls  1.27  cm  thick.  These  two  ducts  were  joined  by  flanges  to  the  bend,  with 
special  care  taken  to  avoid  possible  mismatches  between  the  component  cross- 
sections  which  otherwise  might  disturb  the  flow. 

A  flow-straightening  section  7.15  hydraulic  diameters  long  was  placed  up¬ 
stream  of  the  straight  duct  attached  to  the  bend  inlet  plane.  The  purpose  of 
this  device  was  to  uniformize  and  accelerate  the  development  of  the  cross- 
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stream  plane  distribution  of  the  flow  approaching  the  bend  in  the  upstream 
tangent.  Various  arrangements  of  differently  sized  stainless  steel  screens 
were  tested  in  combination  with  one  or  more  perforated  plexiglass  plates 
3.175  cm  thick,  containing  85  holes  of  3.175  mm  diameter  arranged  in  a 
rectangular  array  spaced  4.495  mm  on  the  centers  in  each  direction.  The 
most  successful  arrangement  of  plates  and  screens  in  the  uniformizing 
section  was  found  experimentally  and  is  shown  in  Figure  2. 

The  test  section  was  part  of  a  closed  loop  system  through  which  water 
at  20°C  was  made  to  flow  by  gravity  from  a  constant  head  tank.  From  this 
tank  the  flow  passed  through  the  test  section,  and  then  into  a  large  sump 
tank  from  where  it  was  pumped  back  to  the  constant  head  tank  by  a  3/4  HP 
Burkes  centrifugal  pump.  The  constant  head  was  ensured  by  a  large  diameter 
PVC  overflow  pipe  extending  through  the  bottom  of  the  head  tank.  Flow  to  the 
head  tank  was  controlled  by  a  gate  valve  and  measured  using  a  2  in.  Barco 
venturi  meter  (P.N.  2-393)  connected  to  a  50  in.  (1.27  m)  differential  mercury 
manometer.  In  order  to  rule  out  the  possibility  of  propagating  perturbations 
induced  by  flow  components,  the  use  of  valves,  sharp  bends  and  metering  devices 
was  avoided  altogether  along  the  test  section  flow  loop.  Flow  to  and  from  the 
test  section  tangents  was  channelled  through  2  in.  i.d.  tygon  tube  pieces;  flexi¬ 
ble  enough  to  be  bent  without  kinks  over  a  large  radius  of  curvature,  yet  stiff 
enough  to  avoid  wall  collapse  due  to  flow- induced  pressure  drop.  Baffles 
located  in  the  constant  head  tank  served  to  dampen  the  swirling  motion  of  the 
flow  leaving  the  tank.  Residual  swirl  from  the  head  tank  and  secondary  motions 
induced  by  mild  curvature  in  the  tygon  tube  upstream  of  the  test  section  were 
eliminated  by  the  flow- uniformizing  section  placed  between  the  tygon  tube  and 
the  upstream  tangent.  All  experiments  were  conducted  for  the  flow  rate  condi¬ 
tion  Imposed  by  the  constant  head  tank.  This  corresponded  to  a  Reynolds  number 
of  Re  *  62,000  and  a  Dean  number  of  De  *  24,000  in  the  flow  test  section. 
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Measurements  of  the  mean  flow  and  turbulence  characteristics  were  made 
using  the  laser-Doppler  velocimeter  technique  in  backscatter  mode  [20 ] -  The 
velocimeter  employed  is  shown  schematically  in  Figure  3  in  relation  to  the 
flow  test  section.  It  comprised  a  2  watt  Lexel  Argon- Ion  water-cooled  laser, 
a  mirror  stage  for  reflecting  the  laser  beam  180°  into  the  velocimeter  optics, 
the  optics,  and  a  4  in  (10.2  cm)  diameter  mirror  for  reflecting  the  converging 
velocimeter  beams  from  the  horizontal  to  the  vertical  direction.  This  mirrror 
also  served  the  function  of  reflecting  the  back-scattered  Doppler-shifted 
radiation  into  the  velocimeter  collecting  optics.  The  optics  were  of  the  DISA 
55X  Modular  Series  and  consisted  of:  two  separately  adjustable  quarter-wave 
retardation  plates;  a  50:50  neutral  beam  splitter;  a  beam  color  splitter;  a 
backscatter  unit  (containing  a  mirror  inclined  45°  with  respect  to  the  veloci¬ 
meter  optical  axis,  and  serving  to  support  at  right  angles  to  the  optical  axis 
the  photomultiplier  optics  consisting  of:  a  color  separator;  two  interference 
filters  and  two  RCA-4526  photomultiplier  tubes);  a  pinhole  section;  a  beam 
translator  and  a  310  mm  achromatic  focusing  lens. 

The  laser  and  velocimeter  optics  were  mounted  to  the  top  of  a  thick 
aluminum  table,  which  was  itself  firmly  bolted  to  an  x,  y,  z  traversing 
mechanism.  The  traversing  mechanism  could  displace  the  table  top  ±  7.5  cm  in 
5  pm  increments  along  any  of  the  coordinate  axes  by  means  of  three  linearly 
encoded  stepping  motors  monitored  by  the  Digital  Equipment  Corporation  PDP  11/34 
minicomputer.  The  minicomputer  functioned  as  the  central  data  acquisition  and 
reduction  controller.  In  addition  to  directing  the  spatial  sequence  of  an 
experimental  run,  the  computer  was  programmed  to  conduct  the  acquisition, 
statistical  processing,  plotting  and  storage  of  Doppler  data  validated  and 
measured  by  a  DISA  55L96  Doppler  signal  processor  or  "counter1'.  The  PDP  11/34 
has  a  256  K  16  bit  random  access  memory  and  is  equipped  with  dual  hard  RL01 
magnetic  disc  drives  (5  Mbytes  each).  The  computer  interacts  by  means  of  an 
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RT  11  software  package  with  various  input-output  devices,  including  a 
Tektronix  4025  graphics  terminal,  a  Decwriter  II  hardcopy  terminal  and  a 
Tektronix  4662  digital  pen  plotter. 

2.2  Experimental  Methodology  and  Measurement  Error 

Prior  to  an  experimental  run,  water  was  allowed  to  flow  through  the  rig 
until  it  was  purged  of  air  bubbles  and  had  attained  a  steady  thermal  state 
corresponding  to  20°C  (±  1°C).  Mass  flow  through  the  test  section  was  con¬ 
trolled  by  setting  the  constant  head  overflow  condition  to  a  mere  trickle 
and  continuously  monitoring  the  pressure  drop  through  the  venturi  meter  con¬ 
nected  to  the  head  tank  feed  line. 

At  any  given  streamwise  measurement  station,  the  velocimeter  support  table 
was  manually  positioned  such  that  the  velocimeter  horizontal  optical  axis  was 
oriented  perpendicular  to  the  test  section  side;  see  Figure  3  for  an  example 
corresponding  to  a  streamwise  location  of  0  *  180°.  Fine  adjustments  to  the  90° 
beam  deflector  mirror  ensured  that  the  velocimeter  vertical  optical  axis  was  per¬ 
pendicular  to  the  test  section  top  surface.  In  combination,  these  adjustments 
produced  a  fringe  pattern  which  was  parallel  (to  within  ±  0.3°  at  any  y,  z 
cross-section  location)  to  the  cross-stream  plane  of  the  test  section.  The 
velocimeter  optical  probe  volume  was  formed  by  the  intersection  of  two  514.5  nm 
(green)  light  beams  with  a  half- angle  in  air  of  4.90°,  for  which  the  volume 
characteristics  were:  a  diameter  of  0.09  mm,  a  length  of  1.1  mm,  and  a  fringe 
spacing  of  3.02  urn  with  about  28  fringes  contained  in  the  probe.  In  reality, 
spatial  filtering  and  threshhold  settings  on  the  counter  reduced  the  dimensions 
of  the  optical  probe.  The  probe  volume  was  positioned  at  the  top  outer-wall 
corner.  Inside  the  test  section,  by  fine  control  of  the  motorized  traversing 
table.  Positioning  in  the  x  and  y  coordinate  directions  was  accurate  to  within 
±  0.05  mm  while  positioning  in  the  z  direction  was  accurate  to  better  than  ±  0.5  mm. 
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With  the  reference  corner  position  at  a  streamwise  location  established, 
the  computer  software  was  activated  which  controls  signal  acquisition  and  data 
processing  on  a  sequentially  scanned  measurement-grid.  After  flow  symmetry  had 
been  established  at  various  streanwise  stations,  the  bulk  of  the  measurements  were 
restricted  to  a  symmetrical  half  of  the  test  section,  on  a  grid  consisting  of 
4  to  5  profiles  at  different  z  locations,  each  containing  29  to  31  equally 
spaced  points  in  the  y  direction.  At  each  point  on  the  measurement-grid  the 
mean  flow  and  turbulence  characteristics  were  statistically  determined  from 
populations  of  5  to  10  samples  consisting  of  1,000  measurements  each.  Each 
measurement  was  required  to  satisfy  the  counter  5/8  validation  comparison  to 
within  a  preset  tolerance  of  3%.  At  every  validation  of  a  Doppler  burst  a 
"data  ready"  signal  was  issued  by  the  counter  to  a  logic  conversion  circuit. 

This  circuit  then  sent  a  triggering  pulse  to  the  computer  parallel  line 
interface  module  which  was  checked  for  data  availability  by  a  software  loop 
approximately  every  20  ys.  No  Interrupt  routines  were  used  for  obtaining 
data  due  to  the  higher  sample  rates  made  possible  by  the  handshake  technique. 

Data  rates  of  about  1  kHz,  with  approximately  60%  validation,  were  obtained 
after  the  flow  was  seeded  with  cornstarch  particles  ranging  In  size  between  1 
and  10  ym.  Optical  alignment  and  automatic  grid  scanning  were  performed  at 
the  following  streamwise  locations:  XH  =  -5,  -1,  1,  5,  10,  20  in  the  straight 
ducts,  and  9  *  3®,  45°,  90®,  130®,  177®  in  the  bend. 

Although  capable  of  two-component  measurements;  the  availability  of  only 
one  counter  during  the  early  stages  of  this  work  restricted  the  use 
of  the  velocimeter  to  single  channel  mode  for  most  of  the  experiments. 

Values  for  the  streamwise  component  of  mean  velocity  and  normal  stress, 

U,  and  u0,  were  derived  directly  from  measurements  obtained  with  the  veloci¬ 
meter  fringes  aligned  perpendicularly  to  the  streamwise  coordinate  direction. 
Values  for  the  transverse  components,  Ur  and  ur,  and  for  the  Reynolds  stress. 
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uQu„,  were  derived  as  described  in  [2,5]  by  combining  measurements  obtained 

with  fringes  oriented  at  angles  of  +  45°  and  -  45°  with  respect  to  the  stream- 

“2 

wise  direction.  As  a  check,  measurements  of  U.  and  u.  obtained  in  this 

0  0 

manner  were  found  to  agree  to  within  experimental  error  with  the  direct 
observations.  Measurements  involving  the  velocity  component  in  the  spanwise 
(z)  coordinate  direction  could  not  be  made  accurately  due  to  optical  inaccessi¬ 
bility  of  the  flow. 

While  the  bulk  of  the  measurements  were  made  in  the  manner  described 
above,  with  the  velocimeter  in  single  channel  mode,  towards  the  end  of  the 
experiment  the  availability  of  a  second  counter  allowed  us  to  perform  a 
limited  number  of  measurements  using  the  two  velocimeter  channels  simultan¬ 
eously.  For  this  case  the  interference  fringe  patterns  were  respectively 
aligned  parallel  and  perpendicular  to  the  streamwise  velocity  component.  To 
resolve  flow  directional  ambiguity,  and  also  to  optimize  the  filter  range  of 
the  counters,  a  net  frequency  shift  of  700  kHz  was  imposed  on  both  channels 
using  a  DISA  55N10  Brag  cell  combined  with  electronic  downmixing.  In  this 
way  additional  measurements  at  *  -5,  -1,  1  and  5  in  the  respective  tangents 
were  obtained.  In  addition,  careful  checks  were  conducted  of  the  data  pre¬ 
viously  obtained  at  =  -5  and  5  and  at  9  =  3°  and  177°.  In  general,  the 
checks  showed  very  good  agreement  between  the  two  measurement  methods,  for 
both  of  the  velocity  components  and  their  respective  normal  stresses.  However, 
at  X^  *  -5  the  two  methods  differed  markedly  with  respect  to  the  measurement 
of  the  very  weak  transverse  velocity  component  arising  at  this  location.  The 
first  method  consistently  showed  larger  scatter  in  the  measurements  and  the 
uncertainty  was  ultimately  traced  to  relatively  small  but  significant  inaccura¬ 
cies  associated  with  determining  the  ±  45°  orientations  required  by  this  method. 
At  =  -5,  -1,  1  and  5,  and  =  3°,  it  is  the  more  accurate  frequency- shifted 
transverse  velocity  data  which  is  reported  here. 


Error  sources  affecting  the  accuracy  (systematic  error)  and  precision 
(random  error)  of  laser-Doppler  measurements  have  been  discussed  by.  for 
example.  Durst,  Melling  and  Whitelaw  [20],  Drain  [21]  and  Buchhave  [22].  In 
this  study  the  most  serious  systematic  errors  were  attributed  to  velocity 
gradient  broadening  and  velocity  bias  respectively.  Velocity  gradient 
broadening  has  been  analyzed  by  Melling  [23]  who  proposed  a  simple  method 
for  estimating  its  magnitude.  Various  weighting  methods  have  been  proposed 
to  correct  for  the  velocity  bias  effect  [21-22],  but  none  of  these  is  entirely 
satisfactory;  they  all  involve  assumptions  regarding  the  statistical  distri¬ 
bution  of  particles  in  the  flow  and,  in  practice,  the  corrections  can  be 
influenced  by  the  additional  problem  of  "incomplete- signal  bias"  [21-22]. 

For  the  conditions  of  this  study,  gradient  broadening  and  velocity  bias  were 
estimated  to  be  significant  only  in  the  near  wall  regions  of  the  flow,  where 
velocity  gradients  and  turbulence  intensities  were  highest  and  the  data  rate 
lowest.  Fortunately,  the  errors  are  of  similar  magnitude  and  of  opposite  sign, 
tending  to  cancel  out  their  respective  effects  on  the  measurements.  For  this 
reason,  and  because  the  errors  were  small  anyway,  corrections  were  not  applied 
to  the  measurements.  Table  1  presents  estimates  of  the  maximum  combined 
inaccuracies  of  these  two  error  sources  on  the  quantities  measured. 

Values  of  the  transverse  mean  velocity  component,  Ur,  and  of  the  Reynolds 

stress,  u.u.  were  prone  to  a  third  systematic  error.  At  the  last  four 
o  r 

locations  in  the  bend,  and  at  XH  *  10  and  20,  these  two  quantities  were 
derived  from  measurements  taken  at  +  45°  and  -  45°  to  the  streamwise  flow 
direction.  An  error  in  setting  this  90°  angle  could  seriously  affect  the 
accuracy  of  these  measurements.  The  problem  has  been  considered  by  Humphrey 
[24]  for  a  flow  of  similar  characteristics  to  the  present  one,  and  he  shows 
that  an  angular  uncertainty  of  0.4°  can  lead  to  an  error  of  5.,  in  u  and  t  3Z 


in  u9ur.  In  this  work  special  care  was  taken  to  ensure  a  maximum  angular 
uncertainty  of  less  than  0.3°  in  setting  the  ±  45°  beam  orientations. 

The  two  main  sources  of  random  error  affecting  the  precision  of  the 
measurements  were  attributed  to  statistical  sampling  uncertainty  (due  to 
the  finite  size  of  sample  populations)  and  uncertainty  in  the  determination 
of  the  reference  or  normalizing  velocity,  Ug.  Estimates  of  the  first  uncer¬ 
tainty  were  derived  form  the  measurements  themselves  and,  for  all  quantities, 
were  found  to  be  less  than  ±  1%  r.m.s.  error.  The  error  in  Ug  was  larger 
(±  2%)  and  arose  principally  from  uncertainties  in  the  construction  of  the 
venturi  meter.  As  a  further  check,  the  bulk  mass  flow  at  each  longitudinal 
station  was  estimated  by  integrating  the  measurements.  This  yielded  a  value 
of  Ug  which  was  within  ±  6*  of  the  venturi  meter  measurement.  Estimates  of 
the  maximum  combined  effects  of  these  two  errors  on  the  measurements  are 
provided  In  Table  1. 

Measurements  of  the  pressure  coefficient,  Cp,  in  the  straight  and  curved 
duct  test  sections,  using  side  wall  pressure  tappings  connected  to  an  inclined 
manometer  bank,  were  also  prone  to  a  random  r.m.s.  error  ranging  from  ±  10%  at 
low  absolute  values  of  Cp  to  ±  5%  at  the  higher  values.  The  uncertainty  in  Cp 
was  due  mainly  to  the  reading  error  associated  with  the  manometer  bank. 
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3.  THE  NUMERICAL  PROCEDURE 

3.1  Governing  Equations  and  Boundary  Conditions 

Three-dimensional,  time  averaged  Navier-Stokes  equations  in  cylindrical 

coordinates,  modeled  according  to  a  two-equation  (k-c)  model  of  turbulence, 

2  2 

are  given  in  [2].  Except  for  the  streamwise  diffusion  terms  (3  /36  ),  exactly 

the  same  modeled  partial  differential  equations  were  used  for  the  partially- 

parabolic  numerical  predictions  of  this  study.  The  k-e  model  presumes  the 

Boussinesq  assumption,  relating  the  stress  and  strain  fields  through  the  notion 

of  an  isotropic  turbulent  viscosity,  v^.  This  quantity  if  determined  from  the 
2 

relation  v.  =  C  k  /e  where,  to  a  first  approximation,  C  is  a  constant  equal 

t  p  P 

to  0.09. 

The  specification  of  boundary  conditions  was  as  follows.  At  the  inlet 

plane  (XH  a  -5)  Ux>  U  and  U2  were  prescribed  from  measurements  of  this  work; 

~7  ~7 

values  of  k  were  estimated  from  measurements  of  u  and  u  (using  symmetry 

y 

considerations  to  derive  uz,  as  was  done  for  Uz);  and  values  for  e  were  speci¬ 
fied  as  in  [25]  by  setting: 

lm  *  <Dh[8x(1  -  1.2  6X)  •  Sy(l  -  1.2  By)]1/2  (1) 


with 


in 


and  S  *-  f  (1  -  rf) 
y  uH  uH 


e 


r  3/4  k3{£ 
w  m 


CQ  =  0.08 
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At  the  exit  plane  (XH  =  5)  two  types  of  boundary  conditions  were  tried.  First 
3P/3x  was  continuously  adjusted  to  conform  with  streamwise  continuity  of  mass, 
while  and  g-(3P/3z)  were  set  e^ua^  t0  zer0-  However,  the  results 

obtained  did  not  differ  from  those  calculated  by  setting  3U/3x  =  0  at  the 
exit  plane  instead.  Because  the  k-e  model  applies  only  to  high  Reynolds 
number  flows,  special  treatment  of  the  flow  near  the  walls  is  required.  This 
is  achieved  by  assuming  the  existence  of  a  logarithmic  velocity  distribution 
between  any  wall  and  the  computed  nodes  nearest  that  wall.  The  additional 
assumption  of  an  equilibrium  flow  in  this  region  leads  to  a  specification 
of  the  wall  shear  stress  in  terms  of  the  wall -node  velocity,  kinetic  energy 
and  logarithmic  law  constants.  Special  treatments  are  also  required  for  near 
wall  determinations  of  k  and  e,  and  are  described  in  [2].  It  should  be  noticed 
that,  except  for  the  manner  in  which  e  was  specified  at  the  inlet  plane  boundary, 
the  turbulence  model  in  this  study  was  identical  to  that  in  [2]. 

3.2  Finite  Differencing  and  Numerical  Algorithm 

Finite  difference  forms  of  the  transport  equations  appropriate  to 
partially-parabol ic  flows  were  derived  following  the  guidelines  given  by 
Pratap  [26].  As  in  [26],  an  upwind  differencing  scheme  was  retained  for  the 
streamwise  velocity  component.  However,  the  second-order  accurate  quadratic 
upstream  interpolation  technique  (QUICK)  proposed  by  Leonard  [27]  and  tested 
by  Han  et  al  [12]  for  turbulent  flow  was  included,  providing  a  choice  for  the 
finite  differencing  of  all  convection  terms  contained  in  the  cross-stream  planes. 

The  QUICK  scheme  was  optional  to  the  k  and  e  equations  as  well  as  to  the  equations 
for  velocity  components.  The  implementation  of  the  scheme  in  the  numerical  algo¬ 
rithm  is  summarized  in  Appendix  I,  a  more  detailed  exposition  is  given  by  Chang  [28]. 

The  elliptic  numerical  procedure  of  [2]  was  modified  to  perform  the 
partial ly-parabol ic  calculations  of  this  study.  The  general  form  of  the 
finite  difference  expression  solved  by  the  modified  procedure  is  given  by: 
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4>p  =  A|  c}>E  +  Aj  4>w  +  Aj  4>n  +  A|  bs  +  B*  (+  D  AP)^  (3) 

where  <j>p  represents  any  one  of  the  dependent  variables  solved  at  a  position 
P  in  the  discretized  flow  domain.  The  coefficients  are  determined  at 
their  respective  (East,  West,  ...)  cell  surfaces,  and  represent  the  combined 
contributions  of  convection  and  diffusion  to  the  balance  of  <*>.  The  term 
contains  contributions  arising  from  sources  (or  sinks)  in  the  flow,  plus  the 
upstream  convective  contribution  to  the  balance  of  <t>.  The  term  in  parentheses 
in  Equation  1  arises  only  when  <p  represents  a  velocity  component,  and  accounts 
for  the  pressure  gradient  contribution  to  the  balance  of  momentum.  As  in  [2], 
the  solution  procedure  employs  a  guess-and-correct  sequence  for  the  calculation 
of  the  three-dimensional  pressure  field.  This  involves  the  iterative  calculation 
of  both  velocities  and  pressure,  which  are  continuously  corrected  until  the 
velocity  components  are  in  conformity  with  local  (and  global)  mass  continuity. 

In  the  present  scheme,  however,  the  flow  variables  are  computed  by 
marching  through  the  flow  domain  along  the  main  flow  direction,  with  the  flow 
variables  at  any  streamwise  plane  being  calculated  solely  from  values  at  the 
upstream  plane.  The  pressure  field  is  always  stored  as  a  complete  three-dimen¬ 
sional  array.  The  remaining  variables  and  coefficients  are  continuously  recal¬ 
culated  and  are  stored  temporarily  as  two-dimensional  arrays  only  at  the  com¬ 
puted  plane  and  the  upstream  plane  adjacent  to  it.  Details  of  the  calculation 
sequence,  including  the  pressure  and  velocity  correction  procedure  and  the  tri- 
diagonal  matrix  algorithm  solution  technique,  will  be  found  in  [26]. 

3. 3  Testing 

The  correctness  of  the  numerical  procedure  and  of  the  turbulence  model 
were  separately  checked  by  reference  to  calculations  of  several  appropriately 
documented  flows.  The  predicted  test  cases  were:  2D,  curved  channel,  developing 
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laminar  and  turbulent  flows  [7];  3D  straight  duct,  developing  laminar  flows 

with  [12]  and  without  [29]  a  Wall  sliding  at  right  angles  to  the  flow;  3D 

straight  duct,  developing  turbulent  flew  [30];  3D,  90°  bend,  developing 

\ 

laminar  flows  with  thin  [3]  and  thick  [24]  inlet  boundary  layers;  and  3D, 

90°  bend,  developing  turbulent  flow  [2].  Except  for  the  last  configuration, 
agreement  was  excellent  between  measurements  and  predictions  of  mean  flow 
quantities  for  all  the  test  cases.  As  anticipated,  the  weak  Reynolds  stress 
driven  secondary  motions  in  the  flow  of  [30]  were  not  predicted  by  the  k-e 
model . 

Due  to  space  limitations,  detailed  comparisons  between  measurements  and 
predictions  of  the  test  cases  are  not  provided  here  but  are  available  in  [28]. 
For  the  case  of  turbulent  flow  in  a  90°  bend,  comparisons  have  been  prepared  by 
Chang  et  al  [11]  for  the  1981-1982  AF0SR-HTTM- Stan ford  Conference  on  Complex 
Turbulent  Flows.  Predictions  of  this  flow  using  the  QUICK  scheme  show  very 
encouraging  agreement  with  corresponding  measurements  over  an  initial  bend 
angle  of  45°,  but  as  of  this  position  discrepancies  of  the  order  of  10-20% 
appear  between  measured  and  predicted  values  of  U0.  In  particular,  the 
transverse  location  of  maximum  UQ  on  the  symmetry  plane  at  9  =  90°  is  placed 
at  n  1  0.82  by  the  calculations,  but  is  located  at  0.60  in  the 
measurements.  Notwithstanding,  these  discrepancies  are  smaller  than  those 
found  in  [2]  using  the  HYBRID  (upwind/central)  differencing  technique  and  an 
elliptic  numerical  procedure. 

Comparisons  between  the  test  cases  predicted  using  both  the  HYBRID  and 
the  QUICK  schemes  showed  clearly  the  superior  performance  of  the  latter  for 
equivalent  grid  sizes  and  distribution  of  nodes.  For  example,  calculations 
of  the  90°  bend  laminar  flow  test  cases  using  the  QUICK  scheme  on  a 
(z  s  11)  *  (r  =  17)  x  (9  =  36)  grid  were  as  accurate  as  corresponding 
calculations  using  the  HYBRID  scheme  on  a  (z  =  15)  ■<  (r  =  25)  •  (fi  =  36) 
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grid.  The  predictions  presented  here  for  turbulent  flow  were  restricted  to 
an  unequally  spaced  grid  consisting  of  the  following  maximum  distribution 
of  grid  nodes:  (z  =  14)  *  (r  =  25)  x  (e  =  45)  in  the  180°  bend,  and 
(z  -  14)  x  (y  =  25)  x  (x  =  20)  in  the  upstream  and  downstream  tangents, 
both  of  length  5  D^.  A  typical  converged  run  time  for  this  grid  represents 

_5 

10  CPUs  per  node  visitation,  and  required  151  kg  words  of  storage  and  a 
total  run  time  of  1600s  for  50  Iterations.  The  criterion  for  convergence  was 
that  the  maximum  of  the  normalized  residual  summations  should  be  less  than 
10‘3. 

Calculations  of  the  flow  were  also  performed  on  coarser  grids.  Compari¬ 
sons  of  these  results  with  the  calculations  provided  here  suggest  that  sub¬ 
stantial  grid  refinements  would  be  required  to  significantly  improve  the 
numerical  accuracy  of  the  present  predictions.  Computation  costs  simply 
prohibited  the  use  of  finer  grids,  or  of  extensive  exploratory  searches  to 
optimize  the  distribution  of  grid  nodes  available. 


4.  RESULTS  AND  DISCUSSION 


Measurements  of  turbulent  flow  in  the  180°  bend,  and  in  its  upstream 
and  downstream  tangents,  are  presented  and  discussed  in  this  section. 

These  are  compared  with  corresponding  predictions  at  selected  streamwise 
locations.  Laminar  flow  predictions  for  the  same  configuration  are  also 
presented  and  discussed.  Even  though  laminar  flow  measurements 
are  lacking  for  comparison,  the  predictions  are  of  value  for:  a)  contrasting 
notable  differences  which  arise  in  the  calculations  from  using  the  HYBRID  and 
the  QUICK  schemes  respectively,  without  the  complicating  factors  associated 
with  turbulence  model  inaccuracies;  b)  helping  to  understand  certain  features 
of  the  turbulent  flow  which  also  arise  in  laminar  flow  regime. 

In  the  figures  corresponding  to  turbulent  flow,  calculated  profiles  are 
plotted  for  2z/DH  =  0  and  0.50  only.  Where  the  HYBRID  and  QUICK  results 
coincide,  a  single  profile  is  shown.  Otherwise,  HYBRID  results  are  plotted 
as  continuous  lines  and  QUICK  as  dashed  lines.  The  reader  should  note  that 
dashed  lines  have  also  been  used  in  some  graphs  for  plotting  best  fits  to 
experimental  data  with  a  larger  than  normal  degree  of  scatter  at  the  2z/DH  =  0. 
and  0.875  spanwise  locations.  Because  photographic  reduction  of  the  figures 
presented  in  this  section  for  discussion  has  reduced  the  accuracy  with  which 
readings  can  be  made  from  them,  full  scale  plots  of  the  turbulent  flow  data 
are  provided  in  Appendix  II.  Laminar  flow  calculations  are  given  in  Appendix  II 

4. 1  Turbulent  Flow  Measurements  and  Calculations 

Measurements  of  the  flow  in  the  upstream  tangent,  taken  at  *  -5  and 
-1  are  shown  in  Figures  6  and  7  respectively.  Comparing  the  data  at  these 
two  stations,  particularly  the  Reynolds  stress  measurements  shown  in  Figure  17, 
suggests  that  the  flow  is  still  developing  at  =  -5,  after  travelling  a 
distance  of  30  0H  from  the  uni formi zing  section.  Figures  5  and  17  provide 


comparisons  of  the  measurements  at  XH  =  -1  with  corresponding  profiles  inter¬ 
polated  from  the  data  obtained  by  Welling  and  Whitelaw  [5]  at  36.8  hydraulic 
diameters  in  a  straight  duct  of  square  cross-section.  The  two  data  sets  are 
in  agreement  to  within  the  experimental  error  of  the  measurements.  Since 
the  measurements  by  Melling  and  Whitelaw  correspond  to  an  essentially 
developed  flow,  the  differences  shown  in  Figure  5  for  the  radial  (transverse) 
velocity  component  are  attributed  mainly  to  the  elliptic  influence  of  the 
bend  on  the  flow  in  the  tangent.  Similar  observations  on  the  influence  of 
a  bend  on  its  upstream  tangent  flow  have  been  noted  by  Humphrey  et  al  [2] 
and  Taylor  et  al  [3]  In  the  same  90°  bend  configuration  but  with  differing 
inlet  conditions.  Taylor  et  al  found  that  in  a  turbulent  flow  with  relatively 
thin  boundary  layers  the  bend  influenced  the  measurements  taken  at  XH  =  -  0.75, 
their  furthest  upstream  location.  Measurements  at  XH  s  -  2.5  taken  by 
Humphrey  et  al  with  essentially  fully  developed  flow  conditions  in  the  up¬ 
stream  tangent  suggested  only  a  weak  elliptic  effect  of  the  bend  on  the  flow 
at  this  location.  However,  these  authors  did  not  measure  the  transverse 
velocity  component;  their  commentary  is  based  on  the  unchanged  appearance 
of  streamwise  velocity  and  turbulence  intensity  contours  between  XH  =  -  11 
and  X^|  =  -  2.5.  By  contrast,  for  the  conditions  of  the  present  experiment, 
the  velocity  data  suggests  that  the  favorable  streamwise  pressure  gradient 
along  the  inner  wall  of  the  bend  (see  Figure  4)  induces  a  mean  transverse 
flow  directed  at  the  inner  wall  which  is  already  noticeable  at  X^  *  -5. 

Although  weak  (Ur/Ug  =  0.02),  it  appears  that  the  transverse  flow  induced 
by  pressure  forces  in  the  upstream  duct  overcomes  the  weaker  stress-driven 
cross-stream  motion  (shown  clearly  in  the  measurements  of  Melling  and  Whitelaw), 
and  results  in  Ur  profiles  whose  shapes,  and  variations  of  shape  with  z 
location,  agree  qualitatively  with  a  simple  superposition  of  the  pressure- 
induced  and  stress-driven  cross-stream  flows. 


Calculations  of  Ur  at  0  =  3°  (and  between  this  location  and  =  -5, 
not  given  here)  always  showed  the  flow  moving  towards  the  convex  wall-side 
of  the  test  section.  Because  the  turbulence  model  is  insensitive  to  stress- 
driven  secondary  motions,  this  result,  and  the  calculated  displacement  of  the 
maximum  in  UQ  towards  the  convex  wall,  clearly  support  the  experimental 
finding  that  elliptic  effects  are  transmitted  from  the  bend  into  the  upstream 
tangent  via  the  pressure  field,  farther  than  has  previously  been  observed. 

As  with  the  Ur  velocity  component,  the  discrepancies  shown  between  calculated 
and  measured  turbulent  stresses  at  6  =  3°  arise  from  the  assumption  of  iso¬ 
tropy  in  the  turbulence  model.  To  predict  more  accurately  the  cross-stream 
motion  and  turbulent  stresses  in  the  upstream  tangent  requires  modeling  of 
Important  anisotropic  flow  characteristics  such  as  near-wall  pressure-strain 
effects.  These  effects  cannot  be  dealt  with  in  the  context  of  the  model 
employed  here. 

The  mean  velocity  and  turbulent  stress  data  taken  at  different  stream- 
wise  locations  In  the  bend  are  plotted  in  Figures  8-12  and  18.  The  pressure 
drop  measured  at  the  side  wall  of  the  bend  is  shown  in  Fig.  4.  In  general, 
the  sense  of  the  flow  up  to  0  *  90°  is  in  agreement  with  earlier  observations 
of,  for  example,  Humphrey  et  al  [2]  and  Taylor  et  al  [3].  Closer  Inspection 
of  the  plots  reveals  additional  interesting  features.  The  measurements  of 
the  pressure  coefficient  Cp  show  the  initially  expected  opposing  pressure 
gradient  at  the  concave  wall  and  the  favorable  gradient  at  the  convex.  In 
contrast  to  the  data  measured  in  [3]  for  a  90°  bend,  the  maximum  and  minimum 
values  of  Cp  do  not  coincide  at  the  same  streamwise  location.  In  this  work 
the  value  of  Cp  maximizes  at  the  concave  wall  at  9  =  45°,  and  attains  its 
minimum  value  at  the  convex  wall  at  9  =  177°.  As  of  9  =  45°,  the  streamwise 
pressure  gradient  is  favorable  throughout  most  of  the  bend.  At  o  -  3"  the 
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streamwise  velocity  profiles  all  shown  their  maximum  values  displaced  towards 
the  inner  radius  wall,  due  to  the  favorable  streamwise  pressure  gradient  there. 
The  radial  component  of  velocity  is  everywhere  directed  towards  the  inner  wall 
except  for  a  small  flow  region  about  the  symmetry  plane.  This  region  marks  the 
inception  of  secondary  flow,  driven  by  the  transverse  pressure  gradient  which 
arises  due  to  lateral  curvature  of  the  main  flow  in  the  bend.  The  secondary  flow 
is  more  intense  at  0  =  45°  as  shown  in  Figure  9.  At  9  =  45°,  the  profile  for  Ur 
near  the  side  wall  (2z/DH  =  0.75)  shows  a  large  negative  velocity,  while  the 
profiles  at  the  other  z  locations  are  large  and  positive.  The  sense  of  motion  is 
from  the  convex  to  the  concave  surface  along  the  symmetry  plane,  and  back  to  the 
convex  surface  along  the  side  walls  of  the  bend.  As  of  0  =  45°,  measurements  of 
Ur  at  all  z  locations  show  this  component  always  positive  (directed  from  the  con¬ 
vex  to  the  concave  surface).  This  means  that  the  cross-stream  plane  return  flow 
adjacent  to  a  flat  wall  in  the  bend  is  confined  to  a  narrow  region  less  than 
D^/8  in  width.  A  simplified  mass  balance  at  0  *  90°  suggests  that  the  radial 
component  of  velocity  in  this  narrow  region  varies  between  0  at  the  wall  and  a 
maximum  of  about  0.40  *  Ug  at  the  peak  location.  This  result  is  in  qualitative 
agreement  with  the  observations  of  Taylor  et  al  [3].  The  highest  return  flow 
measured  by  them  (0.40  *  Ug)  was  found  at  2z/D^  =  0.95  and  8  *  60°  in  a  90°  bend 
with  Re  *  40,000  and  Rc/D  =  2.3.  Calculations  of  the  present  flow  gave  a  peak 
value  of  Ur/Ug  a  0.30  at  2z/DH  a  0.95,  at  the  90°  plane. 

Between  9=3°  and  90°,  the  Influence  of  destabilizing  curvature  raises 
the  levels  of  all  the  measured  Reynolds  stresses  at  the  concave  wall  of  the 
bend.  The  effect  is  particularly  noticeable  in  the  plots  for  the  stresses 
at  0  *  45°,  and  has  decreased  by  the  time  the  flow  reaches  the  0  =  90°  plane. 
Similar  observations  have  been  made  by  Esklnazl  and  Yeh  [31]  In  a  two-dimen¬ 
sional  curved  channel  flow,  and  by  Humphrey  et  al  [2]  and  Taylor  et  al  [3]  in 
their  respective  curved  duct  flows.  In  particular,  Eskinazi  and  Yeh 
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demonstrate  the  strong  generation  of  ur  near  a  concave  wall  and  its  correspond- 
suppression  near  a  convex  wall.  While  at  9  *  45°  the  present  flow  is 
already  three-dimensional,  the  results  for  ur/Ug  are  in  qualitative  agree¬ 
ment  with  the  observations  of  Eskinazi  and  Yeh.  Destabilizing  and  stabiliz¬ 
ing  curvature  effects  at  the  respective  concave  and  convex  walls  of  the 

present  flow  are  responsible  for  producing  large  levels  of  anisotropy. 
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For  example,  at  9  *  45°,  u  0/u  „  *  2  near  the  concave  wall. 

u  r 

Wall-dampening  of  radial  fluctuations  in  the  flow,  coupled  with  pressure 

redistribution  and  turbulent  diffusion  [31]  of  energy  between  the  normal 
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stresses,  account  for  the  larger  levels  of  u  0  . 

At  9  *  90°  and  130%  plots  of  the  two  velocity  components  and  their 
corresponding  normal  stresses  show  striking  variations  in  the  radial  coordinate 
direction.  These  take  the  form  of  relatively  large  decelerations  in  the  mean 
flow  components  at  about  n  *  0.4  and  are  accompanied  by  relatively  large 
increases  in  uQ  and  decreases  in  uf  respectively  at  the  same  locations. 

The  streamwise  deceleration  of  UQ  and  Ur  is  due  to  the  "pumping”  of  low  speed 
fluid  from  the  peripheral  region  of  the  duct  into  the  core  of  the  flow. 

Whereas  viscous  effects  provide  the  mechanism  for  flow  retardation  at  the 
walls.  It  is  the  inviscid  mechanism  of  lateral  flow  curvature  which  drives 
the  cross-stream  motion.  The  phenomenom  has  been  analyzed  and  described  in 
depth  by  Hawthorne  [32],  and  Horlock  and  Lakshminarayana  [33]  discuss  it 
extensively  in  the  context  of  turbomachinery  aerodynamics.  In  particular, 
the  latter  authors  provide  general  expressions  for  the  generation  of  stream- 
wise  vorticity  including  the  influence  of  Bernouilli  surface  rotation  and 
viscous  effects.  Hawthorne  shows  that  in  flows  where  the  angle,  <j>,  between 
the  direction  of  normal  acceleration  and  the  normal  to  the  Bernouilli  surface 
(a  surface  of  constant  total  pressure)  is  $  *  0,  streamwise  vorticity  is  induced 
For  a  small  bend  deflection  angle  &  in  a  flow  with  initial  uniform  vorticity 


in  the  plane  of  the  bend,  the  expression  derived  by  Hawthorne  [32]  is: 

5e  *  -  2  n0  e  <« 

The  same  result  was  obtained  by  Squire  and  Winter  [34]  using  a  different 
analytical  approach  and  is  valid  for  Bernouilli  surfaces  which  remain 
undisturbed  during  passage  through  bend  ($  =  constant). 

It  is  important  to  realize  that  the  inviscid  mechanism  acts  over  the 
entire  cross-section  of  the  flow  (provided  $  *  0).  For  the  kind  of  flow  of 
interest  here,  the  distortion  of  Bernouilli  surfaces  during  passage  of  the 
flow  around  the  bend  cannot  be  neglected.  A  first  approximation  accounting 
for  this  distortion  leads  to  an  equation  for  a,  the  angle  by  which  the 
Bernouilli  surfaces  turn.  The  result  derived  by  Hawthorne  [32]  is: 

This  equation  is  analogous  to  that  governing  the  motion  of  a  pendulum  and, 
as  shown  by  Hawthorne,  predicts  that  the  flow  through  a  bend  should  oscillate 
between  a  =  0  and  a  =  tt  with  a  period  (or  bend  angle)  for  a  complete  oscillation 
approximately  equal  to  2ir //  Rc'/t)H .  The  secondary  motion  is  analogous  to  the 
kinetic  energy  of  the  pendulum  and  also  oscillates  with  the  same  period, 
passing  through  zero  after  each  tt/Dh/Rc  radians  of  turn.  For  the  present 
flow  a  complete  oscillation  corresponds  to  about  197°.  This  means  that  at 
about  98°  the  inviscid  mechanism  starts  working  to  oppose  the  original  sense 
of  the  cross-stream  motion,  with  a  maximum  negative  amplitude  at  0  =  148°. 

The  result,  as  shown  in  Figures  11  and  12  for  0  =  130°  and  177°. 
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For  9  =  90°  and  130*,  plots  of  the  Reynolds  stresses  show  large  changes 

at  the  radial  positions  where  and  Ur  have  been  decelerated.  Transport 
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equations  for  the  turbulent  stresses  can  be  obtained  [35],  and  for  u„  and  u„ 

U  i 

they  take  the  forms  given  below: 
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In  the  above  equations,  P*  and  P*  represent  pressure  strain  redistribution 
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terms,  D.  and  denote  the  effects  of  turbulent  diffusion  and  dissipation, 
and  viscous  diffusion  is  neglected.  The  terms  written  out  explicitly  on  the 
right-hand  side  of  these  two  equations  represent  generation  of  the  stress 
component  and,  hence,  of  the  kinetic  energy  of  turbulence,  k.  Analysis  of 
the  normal  stress  equations,  including  the  effects  due  to  pressure  strain, 
turbulent  diffusion  and  dissipation,  while  desirable,  is  hampered  by  the  un¬ 
availability  of  appropriate  experimental  data;  hence,  the  contributions  of 
the  latter  terms  to  the  respective  balances  cannot  presently  be  established. 
Nevertheless,  simplified  analysis  of  the  generation  terms  alone  provides  a 
basis  for  checking  consistency  in  the  measurements  and  can  shed  light  on  the 
behavior  of  the  flow.  In  the  vicinity  of  the  bend  symmetry  plane,  symmetry 
considerations  support  the  notion  that  the  generation  terms  in  Eqns.  (6  and  7) 
Involving  the  z  coordinate  direction  should  be  small  relative  to  those 

involving  variations  in  the  r  and  0  directions.  Approximate  balances  for 
~2  1 

ufi  and  ur  are  then  given  by: 
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(8) 


(9) 


At  90°  uQur  is  positive  everywhere  and,  except  for  a  small  region 

0.2  <  n  <  0.4  where  3UQ/3r  is  negative,  term  A  represents  a  negative 

"o 

contribution  to  the  balance  of  ua.  From  9  *  45°  to  90°  the  streamwise 

velocity  component  UQ  is  decelerated  between  n  -  0.2  and  0.5  by  the  inviscid 

mechanism  explained  above.  The  result  is  for  term  B  in  Eq.  (8)  to  contribute 

positively  (together  with  term  A)  to  the  balance  of  uQ.  At  all  radial 

locations  near  the  symmetry  plane,  term  C  represents  a  reduction  in  the 
*7 

magnitude  of  uQ,  but  because  of  the  small  values  of  the  ratio  Ur/r,  the 

contribution  of  this  term  relative  to  A  and  B  is  small.  As  shown  in  Figure 

(10),  the  net  result  is  to  produce  a  local  increase  in  the  magnitude  of  uQ 

between  n  s  0.4  and  0.5  approximately. 

In  a  similar  manner,  the  localized  minimum  in  Qr  atn=  0.4  in  the  90° 

plane  can  be  explained.  At  all  radial  locations  term  E  in  Eq.  (9)  contributes 

“7 

positively  to  the  balance  of  ur>  although  weakly  near  the  position  of  the 

minimum  In  ur  because  of  the  small  magnitude  of  u0ur  there.  Between  9  =  45° 

and  90°  the  plots  for  Ur  show  that,  as  for  Ug,  Ur  has  also  been  decelerated 

T 

leading  to  a  positive  contribution  of  term  F  to  the  balance  of  uf.  However, 
the  relatively  strong  streamwise  deceleration  in  Ur  Induces  large  radial 
variations  in  its  profiles  producing  regions  in  the  flow  where  '''Ur/?r  is 
large  and  positive.  As  shown  in  Fig.  (10),  this  occurs  between  n  =  0.2  (at 
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2z/Dh  *  0.75)  and  0.7  (at  2z/C-^  =  0)  and  leads  to  a  negative  contribution  of 

“7 

term  D  to  the  balance  of  ur>  Since  term  D  is  the  only  one  contributing 

negatively  to  up  it  must  be  the  cause  for  the  local  minima  in  the  profiles 
*7 

for  ur  at  9  *  90°.  Therefore,  the  effect  of  this  term  in  the  balance  must 

be  large.  It  appears  that  the  two  terms,  B  in  Eq.  (8)  and  D  in  Eq.  (9),  are 

~2  ~2 

the  principle  source  and  sink  for  ug  and  uf  respectively  in  this  simplified 
analysis.  The  gradual  disappearance  of  the  maxima  in  uQ  and  the  minima  in 
ur  between  9  =  130°  and  177°  is  due  to  the  inviscid  oscillatory  nature  of 
the  flow  since,  by  restoring  high  speed  fluid  to  the  core  of  the  flow  (see 
Figures  11  and  12),  the  respective  roles  of  terms  B  and  D  in  Eqns.  (8  and  9) 
are  reversed. 

~7  “7 

The  anisotropic  effects  in  u0  and  ur  discussed  above  can  not  be  repro¬ 
duced  by  the  calculations.  Nor  is  it  surprising  that  the  minima  observed  in 
the  experimental  profiles  for  the  mean  velocity  components  at  e  *  90°  and 
130°  are  not  predicted.  The  assumption  of  local  isotropy  in  the  turbulence 
model  produces  such  large  levels  of  false  physical  diffusion  as  to  preclude 
an  accurate  spatial  resolution  of  the  flow.  For  example,  at  9  *  90°  in  the 
symmetry  plane  of  the  bend,  estimates  of  v^/v  from  the  measurements  suggest 
that  the  ratio  is  less  than  10  between  n  s  0.1  and  0.9.  However,  calculations 
of  this  ratio  at  the  same  z  =  0  location  vary  from  40  (0  <  n  <  0.4)  to  300 
(0.4  <  n  <  0.95).  Although  numerical  diffusion  compounds  the  problem,  the 
fact  that  it  is  the  more  accurate  QUICK  scheme  results  which  show  worse  agree¬ 
ment  with  the  measurements  at  this  location,  strongly  supports  the  contention 
that  the  source  of  false  diffusion  is  more  physical  in  nature  (due  to  the 
model)  than  it  is  numerical  (due  to  the  differencing  scheme).  In  fact,  the 
higher  levels  of  numerical  diffusion  in  the  HYBRID  scheme  distort  physical 
diffusion  in  the  turbulence  model  and  lead  to  the  incorrect  impression 
of  better  predictions  at  various  streamwise  locations! 


29. 

In  general,  the  QUICK  scheme  calculations  of  mean  radial  velocity  show 
better  qualitative  agreement  with  the  measurements  than  the  HYBRID  reuslts. 

Both  schemes  yield  calculated  values  of  k  which  are  in  reasonably  good  agree¬ 
ment  with  estimates  of  k  from  the  measurements;  with  k  approximated  as 

meas 
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kmeas  =  +  ^  ur)-  However>  agreement  is  misleading,  When  summed, 

~2 

the  experimentally  determined  maxima  in  uQ  and  the  minima  in  uf  compensate 

to  yield  fairly  uniform  radial  distributions  of  k  .  To  indicate  clearly 

meas 

the  inadequacy  of  the  model,  the  plots  provided  show  comparisons  between 

—k  — «  £ 

predicted  values  of  ug  and  ur  (calculated  assuming  uq  =  ur  a  3  ^Calc^  anc* 
direct  measurements  of  these  stresses  at  2z/D^  =  0  and  0.5. 

In  several  of  the  plots,  both  of  the  profiles  given  for  UQ  at  2z/DH  =  0 
and  0.5  respectively  show  larger  values  of  this  component  than  were  actually 
measured,  suggesting  different  values  for  the  experimental  and  computed  mass 
flows  in  the  streamwise  direction.  In  fact,  predictions  of  UQ  nearer  the 
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wall,  not  shown  here,  yielded  values  smaller  than  those  measured  and  gave 
the  required  mass  balance. 

Numerical  experiments  showed  that  a  smaller  length  scale  specification 
in  Eq.  (2)  Um  =  0.01  D^)  at  the  inlet  plane  (X^  =  -5)  led  to  considerably 
reduced  levels  of  the  predicted  turbulent  kinetic  energy.  As  a  consequence, 
even  though  there  was  much  worse  agreement  between  predicted  and  estimated 
values  of  k,  in  the  vicinity  of  0  *  90°  profiles  of  UQ  (but  not  of  Ur) 
showed  the  local  minima  displayed  by  the  measurements.  Again,  this  type  of 
partial  agreement  is  misleading  since  it  is  artifically  reduced  physical 
diffusion  (through  the  specification  of  a  larger  dissipation  at  the  .nlet 
plane)  which  yields  the  result.  In  general,  much  better  overall  predictions 
were  obtained  by  prescribing  an  inlet  plane  length  sca»le  variation  as  out¬ 
lined  in  section  3.1 . 


At  0  =  177°  the  flow  in  the  bend  shows  all  the  maxima  in  the  Ug  profiles 
displaced  towards  the  concave  wall  in  the  presence  of  relatively  large  levels 
of  transverse  flow.  Even  though  the  radial  velocity  component  has  been  de¬ 
creasing  steadily  from  about  9  =  90°  due  to  the  oscillatory  inviscid  flow 
mechanism,  there  is  no  evidence  in  the  Ur  profiles  at  any  location  in  the 
bend  of  a  Dean  instability  in  the  vicinity  of  the  concave  wall.  The  instability 
is  of  the  type  found  by  Taylor  in  concentric  cylinder  flows  [36]  and  has  been 
observed  experimentally  by,  for  example,  Cheng  et  al  [37].  It  has  been  pre¬ 
dicted  numerically  only  in  laminar  flow  regime  [38-40].  Conceivably, 
in  a  bend  of  larger  deflection  angle  than  the  present  one,  the  Inviscid  flow 
oscillation  could  also  induce  an  extra  pair  of  vortices  in  the  vicinity  of  the 
concave  wall.  There  is  supporting  evidence  for  this  conjecture  in  Figure  8 
(Station  12)  of  Hawthorne's  paper  [32]  and  in  the  measured  results  of  this 
work  at  XH  *  1.  At  this  location  the  inviscid  flow  mechanism  has  not  ceased 
to  apply  completely  because  of  the  persistence  of  streamline  curvature  in 
the  flow.  Such  a  mechanism  for  generating  an  extra  pair  of  vortices  in  the 
cross-stream  plane  is  distinct  from  the  Dean  instability  which  is  associated 
with  destabilizing  curvature  effects.  At  0  =  177°  both  of  the  normal  stresses 
display  a  surprising  degree  of  uniformity  and  the  shear  stress  is  small, 
implying  that  the  flow  is  well-mixed  and  relatively  isotropic  at  this  location. 
Except  for  Ur,  at  this  station  the  calculations  are  in  very  good  agreement 
with  the  measurements.  Although  inaccurate  in  terms  of  absolute  values, 

QUICK  scheme  predictions  of  Ur  are  vastly  superior  to  those  obtained  with 
the  HYBRID  scheme.  The  QUICK  scheme  faithfully  reproduces  the  measured  pro¬ 
file  curvatures  while  the  HYBRID  scheme  predicts  negative  values  where  none 
exist  in  the  measurements! 
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At  =  1  (Figure  13)  the  secondary  flow  emerging  from  the  bend  has 
experienced  a  drastic  change  in  both  its  magnitude  and  sense.  The  radial 
velocity  component  has  been  reduced  to  less  than  about  4"  and  a  transverse 
flow  appears  near  the  concave  wall,  opposite  in  sense  to  that  existing  within 
the  bend  at  the  same  radial  location.  A  similar  observation  was  made  by 
Taylor  et  al  [3]  in  their  laminar  flow  measurements  at  XH  *  2.5.  As  mentioned 
above,  the  phenomenon  is  attributed  to  the  persistence  of  the  inviscid  flow 
mechanism  downstream  of  the  bend.  Because  of  the  fairly  uniform  and  similar 
levels  of  the  normal  stress  components  measured  at  9  =  177°  and  XH  =  1 ,  it 
is  unlikely  that  gradients  of  these  stresses  are  altering  significantly  the 
intensity  and  pattern  of  the  cross-stream  flow  at  these  stations. 

As  of  XH  =  0  the  flow  in  the  downstream  tangent  reverts  slowly  to  a 
straight  duct  flow  condition,  with  the  maximum  in  the  UQ  component  tending 
towards  the  duct  center.  However,  the  maximum  has  moved  only  a  small  amount, 
from  n  =  0.9  to  a  value  of  about  0.65  over  a  duct  length  of  20  D^.  Similar 
indications  of  developing  straight  duct  flow  also  appear  in  the  profiles 
for  the  longitudinal  turbulence  intensity  O0/UB  and  the  shear  stress  at 
XH  =  5,  10  and  20  (see  Figures  14-16  and  19).  By  contrast,  the  transverse 
components  of  mean  velocity  and  turbulence  intensity  do  not  show  the  same 
rates  of  adjustment  to  the  new  flow  conditions.  The  profiles  for  Ur  indicate 
that  as  of  XH  =  5  a  cross-stream  flow  pattern,  similar  but  much  weaker  to 
that  in  the  bend,  arises  in  the  downstream  tangent.  At  XH  =  2.5,  their 
furthest  downstream  station,  Taylor  et  al  observed  the  same  result  in  their 
turbulent  flow  measurements. 

The  resurgence  of  a  bend-like  secondary  motion  in  the  downstream 
tangent  is  surprising  and,  at  present,  not  fully  understood.  Calculations 
of  the  flow  between  X^  =  0  and  5,  using  the  QUICK  scheme  are  in  qualitative 


agreement  with  the  measurements.  In  particular,  the  transverse  variations 

of  U„  at  2z/Du  =  0  and  0^.5  are  correctly  simulated  by  QUICK  but  incorrectly 
r  n 

simulated  by  HYBRID.  The  persistence  of  a  peak  in  Ug  near  the  convex-wall 
side  of  the  duct  at  2z/DH  =  0.50  (see  Fig.  14)  is  due  to  turbulence  model 
deficiencies .  Because  of  smearing  by  numerical  diffusion,  this  is  not 
reflected  in  the  HYBRID  calculation. 

Predicted  cross-stream  vector  plots  and  streamwise  velocity  contours  of 
the  flow  at  9  =  177°  using  the  two  differencing  schemes  are  given  in  Fig.  20. 
The  results  show  that  the  HYBRID  scheme  predicts  two  large  vortices  per  half- 
symmetry  plane  in  the  bend  while  the  QUICK  scheme  predicts  essentially  one, 
with  a  second,  much  smaller  counter-rotating  vortex,  at  the  convex  wall. 

The  HYBRID  results  are  in  blatant  disagreement  with  the  measurements,  the 
QUICK  results  are  in  qualitative  agreement  and  this  situation  persists  up 
to  the  exit  plane  at  XH  =  5. 
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4.2  Laminar  Flow  Calculations 

Laminar  flow  calculations  were  performed  for  the  same  configuration 
investigated  in  turbulent  flow  and  are  presented  in  Appendix  III.  Two 
cases  were  predicted  corresponding  to  a  Re  =  890  and  De  =  344.  In  the 
first  case  a  plug  flow  inlet  condition  was  specified  for  U  ,  while  in  the 
second  case  a  parabolic  distribution  was  used.  In  both  cases  the  inlet 
plane  cross-stream  velocity  components  were  set  equal  to  zero.  The  other 
boundary  conditions  used  were:  a)  no  slip  and  impermeable  surfaces  at  all 
walls;  and  b)  fully  developed  flow  condition  at  XH  =  5.  The  inlet  plane  was 
fixed  at  XH  =  -5. 

In  the  plots  cross-stream  and  streamwise  velocities  have  been  normalized 
by  the  bulk  average  velocity,  Ug  *  0.02  m/s.  The  outer  (concave)  wall  position 
is  located  at  the  top  of  each  plot,  while  the  inner  (convex)  wall  position  is 
shown  at  the  bottom.  In  each  figure,  the  top  two  plots  correspond  to  QUICK 
scheme  calculations  and  the  bottom  two  correspond  to  HYBRID  scheme  calcula¬ 
tions.  Because  of  flow  symmetry  only  one  half  of  the  test  section  cross- 
section  is  shown  in  the  figures. 

Plug  Flow  Inlet  Condition 

QUICK  scheme  (QS)  and  HYBRID  scheme  (HS)  calculations  of  this  case  show 
clearly  at  0  =  3°  that  the  flow  entering  the  bend  is  influenced  by  two  effects. 
One  is  the  tendency  of  the  streamwise  velocity  component  to  slow  down  along 
the  periphery  of  the  duct  while  it  speeds  up  in  the  center;  the  other  is  the 
inception  of  the  pressure  driven  secondary  motion,  that  is,  generation  of 
streamwise  vorticity  at  the  side  walls  of  the  bend  by  deflection  of  transverse 
vorticity  in  the  streamwise  direction.  Both  schemes  show  higher  values  of  the 
IJ .  component  at  the  convex  wall. 
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At  3  =  45°  the  HS  shows  the  maximum  in  IL  displaced  towards  the  concave 
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wall,  while  the  QS  places  the  maximum  at  the  center  of  the  duct  and  predicts 
slightly  lower  levels  of  secondary  motion.  Larger  levels  of  numerical  diffu¬ 
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sion  in  the  HS  produce  early  thickening  of  the  side-wall  boundary  layers 
necessary  for  inducing  the  secondary  motion. 

At  e  =  90°  both  schemes  show  the  appearance  of  a  second  but  weaker 
vortex  near  the  convex  wall,  counter-rotating  to  the  main  secondary  motion. 

This  vortex  is  attributed  to  a  less  favorable  streamwise  pressure  gradient  at 
the  convex  wall  as  of  8  =  90°.  At  0  =  130°  the  QS  shows  the  convex  wall  vortex 
extended  over  a  larger  region  on  the  convex  wall-side  of  the  flow,  and  at  177° 
shows  it  has  disappeared.  By  contrast,  the  HS  shows  that  the  convex  wall 
vortex  has  disappeared  by  8  =  130°,  and  instead  predicts  a  Dean-type  instabil¬ 
ity  in  the  form  of  a  new  vortex  at  the  destabilizing  concave  wall.  The  concave- 
wall  vortex  predicted  by  the  HS  persists  at  all  streamwise  locations  as  of 
8  =  130°,  but  never  appears  in  the  QS  calculations.  This  may  seem  puzzling  at 
first,  since  concave  wall  vortex  pairs  are  known  to  arise  in  curved  duct  flow 
[37]  and  are  responsible  for  "pumping"  low  speed  fluid  from  the  concave  wall 
into  the  flow  along  the  symmetry  plane.  However,  they  have  been  predicted 
only  for  fully  developed  laminar  flow  conditions  [38-40],  or  in  developing 
laminar  flows  with  initially  thick  side-wall  boundary  layers  [13]  (flows  with 
relatively  large  viscous  effects).  The  larger  levels  of  numerical  diffusion 
associated  with  the  HS  produce  early  thickening  of  the  side  wall  boundary 
layers  and  yield  conditions  favorable  to  the  formation  of  the  concave  wall 
vortex  pair.  Overall,  the  HS  predicts  a  flow  tending  to  full  development 
more  quickly  than  the  flow  predicted  by  the  QS. 


In  both  cases  the  QS  and  HS  predictions  show  drastic  reductions  in  the 
magnitude  of  the  cross-stream  motion  of  the  flow  leaving  the  bend  which  drops 


from  about  0.20  x  Ug  at  e  =  177°  to  about  0.05  x  Ug  at  =  1 .  There 
Is  a  striking  difference  in  the  sense  of  the  transverse  flow  along  the 
symmetry  plane  at  these  two  locations.  Whereas  the  QS  always  shows  this 
flow  moving  radially  outwards,  the  HS  predicts  a  motion  radially  inwards  in 
the  concave  wall  region  of  the  flow. 

Parabolic  Flow  Inlet  Condition 

Specification  of  a  parabolic  inlet  flow  condition  immediately  imposes 
the  presence  of  transverse  vorticity  throughout  the  entire  flow  entering  the 
bend.  This  is  in  contrast  to  the  plug  flow  case  discussed  above,  for  which 
transverse  vorticity  must  be  diffused  and  convected  from  the  side  walls  of 
the  bend  as  the  flow  evolves.  In  this  sense,  the  present  case  has  more  in 
common  with  the  measured  turbulent  flow  than  the  plug  flow  case. 

At  0  =  3°  both  the  QS  and  the  HS  calculations  show  a  weak  radial  velocity 
component  directed  from  the  concave  to  the  convex  wall,  in  response  to  the 
favorable  3P/30  at  the  convex  wall.  However,  the  U0  contours  show  little 
influence  of  the  transverse  flow.  At  0  =  45°  a  strong  secondary  motion  has 
evolved  and  at  the  side  wall  Uf  =  0.65  x  Ug.  The  higher  secondary  motion 
predicted  relative  to  the  plug  flow  case  is  due  to  the  thicker  side-wall 
boundary  layer  in  the  present  flow. 

At  0  =  45°  and  90°  the  QS  predictions  reveal  the  presence  of  a  corner 
vortex  on  the  convex  wall  due  to  viscous  effects.  This  vortex  is  associate 
with  local  separation  and  produces  the  "wall  stall"  observed  on  the  suction 
surfaces  In  cascade  flows  [33].  The  HS  predictions  fail  to  reveal  Its 
presence. 

Between  0  =  90°  and  177°  transverse  variations  of  UQ  predicted  by  the  QS 
are  in  very  good  qualitative  agreement  with  the  turbulent  flow  measurements; 
corresponding  predictions  using  the  HS  are  not.  The  role  of  the  inviscid 


flow  mechanism  is  clearly  shown:  initially  low  speed  flow  from  the  convex 
wall  is  pumped  into  the  core  of  the  flow  (6  »  90°),  the  trend  being  reversed 
subsequently  to  restore  high  speed  fluid  in  the  core  (e  =  130°).  The  oscil¬ 
lating  character  of  the  flow  in  the  bend  is  also  supported  by  QS  predictions 
of  the  cross-stream  velocities.  The  vector  plots  show  that  the  magnitude  of 
the  secondary  motion  is  largest  at  e  s  45°  and  smallest  at  e  =  130°  (in  agree¬ 
ment  with  the  inviscid  flow  analysis  of  Hawthorne  [32]).  In  the  present  case, 
viscosity  dampens  the  oscillation  and  reduces  the  magnitude  of  the  secondary 
motion  in  passage  through  the  bend.  In  a  long  enough  bend,  a  non-oscillating 
asymptotic  flow  condition  should  eventually  be  reached. 

As  In  the  plug  flow  case,  between  e  =  45°  and  130°  the  QS  predictions 
reveal  a  convex-wall  vortex  adjacent  to  the  symmetry  plane.  By  8  *  177°  the 
vortex  has  disappeared.  The  HS  reveals  the  presence  of  this  vortex  between 
0  =  45°  and  130°  only.  Likewise,  the  presence  of  a  Dean- type  vortex  structure 
at  the  concave  wall  Is  revealed  by  both  calculation  schemes  as  of  6  =  130°. 

At  0  ■  177°  the  size  and  intensity  of  the  Dean-type  vortex  predicted  by  the 
QS  are  larger  than  the  values  predicted  by  the  HS.  Numerical  diffusion  in 
the  HS  predictions  is  responsible  for  smearing  out  this  vortex  structure  and 
reducing  its  intensity  In  passage  through  the  bend.  The  combined  pumping  of 
low  speed  fluid  from  both  curved  walls  into  the  core  flow  along  the  symmetry 
plane  has  the  effect  of  Isolating  a  region  of  high  speed  streamwlse  flow 
between  the  side  wall  and  symmetry  plane  of  the  bend. 

Downstream  of  the  bend,  plots  of  ihe  cross-stream  flow  predicted  using 
the  QS  show  the  quick  decay  of  two  cross-stream  vortical  structures  whose 
origins  were  the  pressure  driven  secondary  motion  and  Dean  instability 
respectively  in  the  bend.  Strictly,  these  two  mechanisms  are  not  Independent 
in  the  present  configuration.  Due  to  the  unit  aspect  ratio  of  the  duct,  there 
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exists  a  synergic  interaction  between  them.  HS  calculations  at  the  same 
locations  (X^  =  1  and  5)  show  a  faster  decay  of  the  structures,  which  is 
attributed  to  numerical  diffusion.  It  is  interesting  to  observe  that  the 
sense  of  the  transverse  velocity  component  predicted  by  the  QS  at  =  1 
is  in  complete  qualitative  agreement  with  the  turbulent  flow  measurement  of 
this  component,  shown  in  Figure  13  of  the  text. 

For  both  the  plug  flow  and  the  parabolic  flow  inlet  distributions,  it 
is  important  to  note  that  the  respective  flows  in  the  bend  did  not  reach  a 
fully  developed  state  and  were  significantly  influenced  by  the  presence  of 
the  downstream  tangent. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  new  measurements  provided  here  for  the  mean  flow  in  strongly 
curved  180°  bends  reveal  features  which  are  in  qualitative  agreement  with 
results  obtained  from  inviscid  flow  analysis.  Measurements  of  the  turbulence 
characteristics  of  the  flow  at  8  =  90°  and  130°  show  striking  variations  in 
the  radial  direction  which  appear  to  be  due  to  large  shearing  motions  induced 
by  the  inviscid  mechanism  in  the  core  of  the  flow.  Within  the  context  of  an 
Interpretation  which  neglects  turbulence  diffusion  and  pressure  redistribution 
effects,  the  velocity  gradients  Induced  in  the  core  of  the  flow  work  locally 
on  the  turbulent  stresses  to  increase  the  longitudinal  stress  component  (uQ) 
and  decrease  the  radial  stress  component  (ur). 

Between  e  =  0°  and  45°,  stabilizing  and  destabilizing  curvature  effects 
at  the  convex  and  concave  wall  respectively  are  responsible  for  generating 
relatively  large  levels  of  anisotropy  near  these  surfaces.  As  of  0  =  45° , 
the  cross-stream  motion  contributes  to  the  production  of  a  more  complicated 
anisotropic  pattern  in  the  turbulent  stresses.  For  example,  regions  of 
negative  production  of  turbulent  kinetic  energy  arise  at  9  =  90°.  Between 
0  =  130°  and  177°  the  anisotropy  In  the  measured  normal  stresses  Is  reduced, 
due  to  the  uni forml zing  Influence  of  the  secondary  motion  In  the  flow. 

Measurements  of  Ur  in  the  bend  show  that,  except  for  around  0  *  45°, 
the  return  flow  In  the  cross-stream  plane  is  restricted  to  fairly  narrow 
bands  of  width  =  D^/8  adjacent  to  each  of  the  side  walls.  The  measurements 
of  Ur  did  not  reveal  vortex  structures  at  the  convex  wall  in  the  bend,  nor 
at  the  concave  wall  where  they  might  be  expected  to  arise  due  to  the  Dean 
Instability.  However,  the  sense  of  the  cross-stream  flow  at  XH  *  1  does 
suggest  a  brief  appearance  of  this  type  of  Instability  in  the  exit  tangent. 
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Between  a  =  0°  and  XH  =  20  there  was  no  evidence  in  the  measurements  of  a 
significant  Reynolds  stress  driven  secondary  motion. 

In  this  investigation  the  influence  of  the  bend  on  the  upstream 
tangent  flow  is  already  noticeable  at  XH  =  -5.  Although  weak,  a  transverse 
flow  is  set  up  at  this  location  (and  at  XH  =  -1)  of  characteristics  which 
correspond  with  the  simple  notion  of  a  superposition  of  two  cross-stream 
motions:  the  first  induced  by  Reynolds  stress  gradients  in  the  cross-stream 
plane;  the  second  due  to  the  favorable  streamwise  pressure  gradient  arising 
at  the  convex  wall  in  the  bend.  In  the  downstream  tangent  secondary  motions 
are  drastically  decreased  between  e  =  180°  and  X^  =  1.  The  decrease  is  re¬ 
lated  to  the  oscillating  nature  of  the  secondary  motion  which,  as  shown  by 
inviscid  flow  theory  and  supported  by  the  present  results,  as  of  e  *  90° 
has  started  opposing  the  Initial  source  of  circulation  set  up  by  the  trans¬ 
verse  pressure  gradient  in  the  bend.  Between  X^  =  1  and  20  the  flow  reestab¬ 
lishes  characteristics  typical  of  straight  ducts.  However,  the  measurements 
show  the  persistence  of  a  weak  secondary  motion  In  the  downstream  tangent 
with  the  same  sense  of  rotation  as  the  flow  in  the  bend  and  a  slow  recovery 
of  the  radial  normal  stress  compared  to  the  longitudinal  component. 

Predictions  of  the  experimental  flow  were  made  with  a  k-e  model  of 
turbulence  using  a  partially-parabolic  numerical  procedure.  Two  finite 
differencing  techniques  were  used  with  global  accuracies  of  first  order 
(HYBRID)  and  second  order  (QUICK).  Testing  and  predictions  with  QUICK  in 
the  laminar  flow  regime  clearly  showed  it  to  be  superior  to  HYBRID.  The 
differences  between  the  results  obtained  using  the  two  schemes  in  turbulent 
flow  are  notable.  Of  the  two,  only  QUICK  reproduces  correctly  the  trends 
In  the  measured  cross-stream  flow.  The  HYBRID  scheme  gives  the  incorrect 
Impression  of  better  streamwise  component  velocity  predictions.  This  is 
due  to  numerical  diffusion  which  smears  out  the  spacial  variation  of  this 
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component.  Although  QUICK  scheme  predictions  for  UQ  differ  more  from  the 

U 

measurements,  this  is  explained  by  the  fact  that  with  numerical  diffusion 
reduced,  it  is  the  false  physical  diffusion  in  the  k-e  model  which  is  re¬ 
vealed  by  the  computed  results. 

Numerical  experiments  confirmed  a  fairly  strong  sensitivity  of  the  model 
to  the  inlet  plane  specification  of  dissipation.  Fixing  a  small-but  incor¬ 
rect-dissipation  length  scale  at  the  inlet  plane  led  to  somewhat  improved 
predictions  of  the  UQ  component  at  0  *  90°,  where  this  component  shows  large 
changes  in  the  transverse  direction.  However,  predictions  of  Ur  and  of  k 
were  much  less  satisfactory. 

In  an  earlier  study  [11]  it  was  reported  that  predictions  of  turbulent 
flow  in  a  90°  bend  were  mildly  sensitive  to  the  sense  of  cross-stream  flow 
fixed  at  the  inlet  boundary  condition.  The  effect  was  too  small  to  alter 
significantly  the  streamwise  evolution  of  the  flow.  In  this  study,  insigni¬ 
ficant  differences  were  observed  between  predictions  of  the  flow  with  the  cross¬ 
stream  motion  set  equal  to  zero,  and  predictions  with  the  cross-stream  motion 
specified  from  the  measurements.  The  insensitivity  of  the  flow  in  a  bend  to 
the  inlet  plane  streamwise  component  of  vorticity  is  explained  by  reference  to 
the  general  form  of  the  equation  for  the  variation  of  this  quantity  along  a 
streamline  [33]: 


In  the  case  of  a  flow  entering  a  bend  from  a  straight  duct,  generally 

£  »  Co  »  especially  near  the  side  walls.  As  a  consequence,  the  equation  shows 
r  w 

that  along  streamlines  curving  through  the  bend,  the  streamwise  development 
of  vorticity  is  due  primarily  to  deflection  of  transverse  vorticity.  Viscous 
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diffusion  of  streamwise  vorticity,  even  in  turbulent  flows,  is  of  secondary 
importance. 

Laminar  flow  calculations  of  the  experimental  configuration  have  been 
performed  for  two  entrance  flow  conditions  corresponding  to  plug  flow  and 
parabolic  flow  profiles.  Both  QUICK  and  HYBRID  scheme  predictions  were 
made.  Significant  differences  between  the  two  schemes  are  attributed  to 
the  larger  levels  of  numerical  diffusion  in  the  HYBRID  differencing  approach. 
QUICK  scheme  calculations  reveal  four  types  of  cross-stream  motion  which 
overlap  at  different  locations  in  the  bend:  the  main  cross  stream  flow, 
driven  by  the  transverse  pressure  gradient;  a  small  vortex  located  at  the 
convex  wall,  of  brief  existence,  associated  with  a  decrease  in  the  favorable 
streamwise  pressure  gradient  at  this  wall  in  passage  through  the  bend;  a 
larger  vortex  located  on  the  concave  wall,  of  a  persistent  nature,  associated 
with  the  Dean  type  instability  arising  in  connection  with  destabilizing 
curvature  effects;  a  small  vortex  located  in  the  convex  wall  corner,  of 
brief  existence,  driven  by  the  side-wall  boundary  layer  and  associated  with 
the  phenomenom  of  stall  in  turbine  cascades.  Oscillation  of  the  main  flow 
was  predicted  numerically  by  the  QUICK  scheme  and  was  most  Intense  in  the 
parabolic  flow  case.  Because  of  the  larger  net  content  of  transverse  vorti¬ 
city  at  the  inlet  plane,  the  parabolic  flow  predictions  resemble  the  turbulent 
flow  measurements  more  than  the  plug  flow  case. 

The  following  recommendations  are  offered  for  continuing  work: 

1.  Laminar  flow  measurements  should  be  performed  in  the  bend  and 

downstream  tangent  to  establish  definitively  the  superior  performance 
of  QUICK  for  predicting  curved  duct  flows  and  for  investigating  the 
various  vortical  structures  predicted  by  this  scheme. 

More  detailed  turbulent  flow  measurements  are  required  to  understand 
the  relaxation  processes  taking  place  in  the  downstream  tangent. 


2. 
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Quanti ty 

Range  of  Maximum 
Systematic  Error  + 
(deviation  %) 

Maximum  ++ 
Random  Error 
(r.m.s.  %) 

VUB 

-  1*  to  2% 

±  2% 

VUB 

-  U  to  2% 

±  2%  to  ±  3% 

56/UB 

-  2%  to  2% 

±  2%  to  ±  3% 

VUB 

-  2%  to  2% 

±  n  to  ±  4% 

VV/UB 

-  2%  to  2% 

±  3%  to  ±  7% 

+  Error  confined  mainly  to  near  wall  flow  regions; 

does  not  include  angular  uncertainty  discussed 
in  text. 

++  Largest  errors  confined  to  small  flow  regions;  does 

not  include  a  small  positioning  uncertainty  of  probe 
volume  in  the  flow. 


Table  1;  Tabulation  of  Estimated  Maximum  Measurement  Errors 


Schematic  showing  test  section  dimensions,  coordinate  system  and 
velocity  components  of  the  flow. 


Figure  2  Optimum  experimental  arrangement  of  stainless  steel  screens  and 
perforated  plate  in  the  flow-uni forml zing  section  preceding  the 
upstream  tangent. 
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Figure  3  Top  and  side  view  of  laser-Doppler  velocimeter  aligned  for 
measurements  at  a  bend  angle  0  *  180°.  Traversing  mechanism 
details  are  not  shown. 
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Figure  4  Distribution  of  pressure  coefficient  (Cp  *  AP/pUg)  In  the  bend 

and  tangents  for  conditions  of  the  experiment:  Re  *  64,000, 

Rc/Dh  *  3.35  and  De  *  24,000.  Measurements  made  through  pressure 

taps  on  side  walls  of  test  section  components. 
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Comparison  between  turbulent  flow  measurements  of  mean  velocity 
and  normal  stress  from  this  work  (o)  at  *  -1  and  measurements 

by  Mel  ling  and  Whitelaw  ( - )  after  a  development  length  of 
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Figure  8  Turbulent  flow  measurements  and  calculations  of  mean  velocity  and 
normal  stress  in  the  bend  at  6  *  3°:  ( — )  HYBRID,  (---)  QUICK; 

dashed  line  at  0.875  is  best  fit  to  measurements. 
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I  Figure  9  Turbulent  flow  measurements  and  calculations  of  mean  velocity  and 

normal  stress  in  the  bend  at  0  *  45°:  ( — )HYBRID,  ( — )  QUICK. 
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Figure  10  Turbulent  flow  measurements  and  calculations  of  mean  velocity  and 

normal  stress  In  the  bend  at  9  =  90°:  ( - )HYBRID,  (  — )  QUICK; 

dashed  line  at  0.875  is  best  fit  to  measurements. 
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Figure  12  Turbulent  flow  measurements  and  calculations  of  mean  velocity  and 

normal  stress  in  the  bend  at  9  =  177°:  ( - )  HYBRID,  (  — )  QUICK; 

dashed  lines  at  0.75  and  0.875  are  best  fits  to  measurements. 
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Figure  14  Turbulent  flow  measurements  and  calculations  of  mean  velocity 
and  normal  stress  in  the  downstream  tangent  at  s  5:  ( — 1 

HYBRID,  ( — )  QUICK;  dashed  line  at  0.875  is  best  fit  to 
measurements. 
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Figure  15  Turbulent  flow  measurements  of  mean  velocity  and  normal  stress 
in  the  downstream  tangent  at  =  10;  dashed  lines  at  0.75  and 

0.875  are  best  fits  to  measurements. 


Figure  17  Turbulent  flow  measurements  of  shear  stress  at  XH  =  -5  and  -1 

(o)  data  from  this  work;  ( - )  data  from  Mel  ling  and  Whltelaw; 

dashed  lines  at  0.75  are  best  fit  to  data. 


r 


(r-  #t> 

XBL  8210-2801 


Figure  18  Turbulent  flow  measurements  of  shear  stress  at  four  longitudinal 
stations  in  the  bend;  dashed  line  at  0.75  is  best  fit  to  measure 
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gure  20  Turbulent  flow  predictions  of  cross-stream  and  streamwise  velocities 
at  8  *  177°:  velocity  components  are  normalized  by  Ug ;  outer 

radius  (concave)  wall  is  shown  at  the  top  of  each  plot. 


APPENDIX  I 


QUICK  SCHEME  FORMULATION  FOR  TRANSPORT  VARIABLES 


1 . 

This  appendix  contains  a  brief  account  of  the  formulation  of  the  QUICK 
scheme  as  well  as  its  incorporation  into  TEACH-type  codes.  It  is  illustrated 
via  a  simple  1-D  problem.  The  appropriate  2-D  counterparts  are  summarized  in 
table  A-l .  Details  of  the  scheme  formulation  can  be  found  in,  for  example, 

[27,28]  and  of  its  testing  in  [12]. 

As  an  example,  consider  the  simplified  convection-diffusion  equation 

&  WW -&(*&)•  0  (A-’> 

Following  the  standard  control  volume  approach  of  TEACH-type  codes  [29], 
equation  (A.l)  is  integrated  over  its  cell,  shown  in  Figure  A.l.  Using 
central  differencing  for  the  diffusion  terms,  equation  (A.l),  with  all  the 
subscripts  and  dimensions  referred  to  Figure  A.l,  becomes 

Ce^e  *  cw*w  "  °e(*E  "  V  +  Dw(l,,P  "  V  “  0  (A-2) 

where  Cg  =  (pU)e 

cw  a  (put. 

°e  5  V{x1 

aid  Dw  V6X1-1 

For  the  convection  terms  C_ <p_  and  C  <j>  ,  Leonard  [27]  proposed  a  quadratic 

6  6  W  W 

upwind  Interpolation  scheme  (QUICK)  to  approximate  4>e  and  $w  at  the  surfaces 
of  the  control  cell.  The  proper  formulae  for  the  east  surface  are: 


*e  s  2^*'E  +  V  '  8'  ^ 


(  '" i )  /^EE  “  *E  *E  ’  *P 


6i+l 


6.  >/A1+l  *  if  Ce<0 


or 


I  ( 6 . ) 

4>e  *  +  M  8  (  6“ 


4>p  - 

-7 - ^  /A, 

Si-1  ’ 


if  C  >  0 
e 


To  facilitate  the  incorporation  of  QUICK  into  equation  (A. 2),  equations 
(A. 3)  are  written  in  the  compact  form: 

<Pe  =  ?(4>E  +  <Pp)  +  ae4>p  +  8e 


where  a  * 1 
e 


'6i(1  +  V6i-1>/8A1 

r  6i/8Ai-i 


if  ce  >  o 


if  Ce  <  0 


and  6g  = 


/-  (6i)Z(4>E/6i  +  V6i-l)/8Ai 


if  ce  >  0 


-  (<5i  )2(4>ee  ”  ^E)/5i+l  ‘  V6i)/8Ai+l  1f  ce  <  0 


Similarly  we  write 


♦w  =  2<*P  *  V  +  Vp  *  8w 


Substituting  equations  (A. 4)  and  (A. 5)  into  equation  (A. 2)  and 


(A. 3) 


rearranging  yields: 


1 


3. 


VDe  ♦  1  Ce  +  -  ?  Cw  -  <\C„  -  “eCe»  ‘  VDe  '  1  Ce> 

+  *W{Dw  *  7  Cw>  *  <Vw  -  6eCe> 

Finally,  adding  the  continuity  relation  Cw  -  Cg  =  0  multiplied  by  $p  to 
both  sides  of  equation  (A. 6),  yields: 

*AE  +  AW  SP^P  =  AE^E  +  +  Su 


(A. 6) 


(A. 7) 


with 


\ 

S„ 


and 


o,  -  7  ce 

D„*K 

a  C  -  a  C 
w  w  e  e 

8  C  -  8  C 
w  w  e  e 


Table  A.l  summarizes  the  appropriate  formulae  for  general  2-D  problems. 

It  is  worth  noting  that  the  separation  of  terms  and  8^  in  Table  A.l  ensures 
that  Sp  is  always  negative.  This  Is  a  highly  desirable  result  for  convergence 
of  the  matrix  inversion  routine. 


Figure  A.l 


TABLE  A.  1 
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APPENDIX  II 


FULL  SCALE  PLOTS  OF  TURBULENT  FLOW  DATA 


Due  to  reproduction  and  mailing  costs,  some  copies  of  this  report  are 
lacking  the  extensive  data  set  pertaining  to  Appendix  II.  If  full-seal 
plots  are  not  included  here,  they  will  be  made  available  upon  request. 
See  title  page  for  mailing  address  of  the  Principal  Investigator. 


gure  1:  Calculated  cross-stream  and  streamwise  velocity  components  at 

=  -1  in  laminar  regime.  Plug  flow  inlet  profile.  Components 

normalized  by  Ug  =  0.02  m/s.  Outer  radius  at  top  of  each  plot. 
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gure  4:  Calculated  cross-stream  and  streamwise  velocity  components  at 

6  =  90°  in  laminar  regime.  Plug  flow  inlet  profile.  Components 
normalized  by  *  0.02  m/s.  Outer  radius  at  top  of  each  plot. 
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gure  5:  Calculated  cross-stream  and  streamwise  velocity  components  at 

9  =  130°  in  laminar  regime.  Plug  flow  inlet  profile.  Components 
normal ized  by  IL  =  0.02  m/s.  Outer  radius  at  top  of  each  plot. 
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Figure  6:  Calculated  cross-stream  and  streamwlse  velocity  components  at 

0  *  177°  in  laminar  regime.  Plug  flow  inlet  profile.  Components 
normalized  by  Ug  =  0.02  m/s.  Outer  radius  at  top  of  each  plot. 
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Figure  7:  Calculated  cross-stream  and  streamwise  velocity  components  at 
XH  =  1  in  laminar  regime.  Plug  flow  inlet  profile.  Components 

normalized  by  UR  =  0.02  m/s.  Outer  radius  at  top  of  each  plot. 


Figure  8:  Calcilated  cross-stream  and  streamwise  velocity  components  at 
XH  =  5  in  laminar  regime.  Plug  flow  inlet  profile.  Components 

normalized  by  =  0.02  m/s.  Outer  radius  at  top  of  each  plot. 
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Figure  9:  Calculated  cross-stream  and  streamwise  velocity  components  at 
0  =  3°  in  laminar  regime.  Parabolic  flow  inlet  profile. 
Components  normalized  by  IL  =  0.02  m/s.  Outer  radius  at  top 
of  each  plot. 
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Figure  10:  Calculated  cross-stream  and  streamwise  velocity  components  at 
9  *  45s  in  laminar  regime.  Parabolic  flow  inlet  profile. 
Components  normalized  by  U„  =  0.02  m/s.  Outer  radius  at  top 
of  each  plot. 
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Figure  11:  Calculated  cross-stream  and  streamwise  velocity  components  at 
9  =  90°  in  laminar  regime.  Parabolic  flow  inlet  profile. 
Components  normalized  by  =  0.02  m/s.  Outer  radius  at  top 
of  each  plot . 
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Figure  14:  Calculated  cross-stream  and  streamwise  velocity  components  at 
XH  ■  1  In  laminar  regime.  Parabolic  flow  Inlet  profile. 

Components  normalized  by  U„  *  0.02  m/s.  Outer  radius  at  top 
of  each  plot.  L 
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Figure  15:  Calculated  cross-stream  and  streamwise  velocity  components  at 
XH  *  5  In  laminar  regime.  Parabolic  flow  Inlet  profile. 

Components  normalized  by  IL  *  0.02  m/s.  Outer  radius  at  top 
of  each  plot. 
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